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ABSTRACT 

We present a study of the infrared properties of X-ray selected, moderate luminosity (i.e, 
Lx= 10 42 — 10 44 ergs s" 1 ) active galactic nuclei (AGNs) up to zss3, in order to explore the 
links between star formation in galaxies and accretion onto their central black holes. We use 
100 \xm and 160 \xm fluxes from GOODS-//e rschel - the deepest survey yet undertaken by 
the Herschel telescope - and show that in the vast majority of cases (i.e., > 94 per cent) these 
fluxes are dominated by emission from the host galaxy. As such, these far-infrared bands 
provide an uncontaminated view of star formation in the AGN host galaxies. We find no 
evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs 
at any redshift, suggesting that star-formation is decoupled from nuclear (i.e., AGN) activity 
in these galaxies. On the other hand, we confirm that the star formation rates of AGN hosts 
increase strongly with redshift; by a factor of from z < 0.1 to z = 2 — 3 for AGNs 

with the same range of X-ray luminosities. This increase is entirely consistent with the factor 
of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (i.e., 
main-sequence) galaxies over the same redshift range. Indeed, the average SSFRs of AGN 
hosts are only marginally (i.e., w 20 per cent) lower than those of main-sequence galaxies 
at all surveyed redshifts, with this small deficit being due to a fraction of AGNs residing in 
quiescent (i.e., low-SSFR) galaxies. We estimate that 79 ± 10 per cent of moderate luminosity 
AGNs are hosted in main-sequence galaxies, 15 ±7 per cent in quiescent galaxies and < 10 
per cent in strongly starbursting galaxies. We derive the fractions of all main sequence galaxies 
at z < 2 that are experiencing a period of moderate nuclear activity, noting that it is strongly 
dependent on galaxy stellar mass (M stars ); rising from just a few per cent at M stars ~ 10 10 M Q 
to >20 per cent at M stars > 10 11 M . Our results indicate that it is galaxy stellar mass that is 
most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue 
that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than 
violent mergers, which suggests that high redshift disk instabilities could be an important 
AGN feeding mechanism. 
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1 INTRODUCTION 

Over the past two decades a significant bank of evidence has de- 
veloped demonstrating, to first order, the close links between the 
growth of galaxies and their resident black holes (e.g., the similar- 
ity between the black hole accretion and star formation histories of 
the Universe, t he tight correlation between black hole and galax 
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Haring & Rix 20QJ; iBovle & Terlevichl 1 19981 ; iFranceschini et al] 
1999; Merloni 2004). Results from some theoretical models sug- 



gest that these links arise through feedback processes between 
the galaxy and its a ccreting black hole, regulating their re- 
spective growths (e. g., iGranato et alj|2004l ; |Pi Matteo et ai] |2005l ; 
IHopkins et al. 2006). However, outside these models the physical 
mechanisms driving this expected feedback are poorly understood, 
primarily due to a lack of constraints from observations. Indeed, 
it is still a matter of debate whether all levels of nuclear activity 
directly impact star formation (or vice-versa), or whether the scale 
of their influence is a function of AGN or galaxy properties (e.g., 
AGN luminosity, galaxy mass, redshift)Q Measuring the simultane- 
ous levels of nuclear and star formation activity taking place within 
the galaxy population is crucial if we are to determine how or, in- 
deed, if these processes are causally linked. 

Measuring the stellar properties of a galaxy (including star 
formation activity) when faced with contamination from an AGN 
presents a significant challenge. The most commonly used diag- 
nostics rely on optical spectral features (e.g., the 4000 A break, the 
strengths of the H84 absorption and [O II] emission lines), which 
have been applied extensively to large scale spectral surveys to 
measure the properties of large numbers of AGN host galaxies (e.g., 
Kauffmann et al. I l2003l ; ISilverman et ai]|2009h . However, spectral 
diagnostics are not immune to AGN contamination and optical di- 
agnostics, in particular, are susceptible to the effects of reddening. 
By contrast, the far-infrared (hereafter, FIR) regime is largely im- 
pervious to obscu ration and, except in the most extreme cases (e.g., 
Hon v et al ■1201 lh . provides a measure of the star formation activity 
of the host galaxy that is largely uncontaminate d by the AGN (e.g., 
iNetzer et al.l2007l ; lHatziminaoglou et al.l201oh . However, it is only 
with the recent availability of Spitzer and, more crucially, Herschel 
that we have been able to measure the FIR emission of AGN hosts 
at high enough redshifts to probe the epochs during which most 
stellar and black hole mass in the Uni verse was built up (i.e., z > 1; 
iMullanev et alj20ld ; lxu"e et al.l20ld ; lshao et al.l201CI) . These first 
FIR studies of AGNs in the high redshift Universe had to contend 
with relatively low infrared detection fractions (typically <20 per 
cent) and, as such, relied heavily on stacking analyses. By contrast, 
in this study we use data from the deepest FIR survey yet under- 
taken by Herschel to reach a FIR detection rate of >40 per cent. 
For the first time, we can use individual detections to explore the 
trends between AGN and star formation activity to high redshifts 
(i.e., z < 3), only resorting to stacking analyses to confirm, rather 
than identify, correlations derived from FIR observations. 

The FIR output of a galaxy provides a direct measure of the 
rate at which that gal axy is forming stars (i.e., the star forma- 
tion rate, or SFR; e.g.. lKennicuttlll998l) . Clearly, SFRs are a cru- 
cial element in most, if not all, studies exploring the links be- 
tween nuclear activity and galaxy growth. However, it has recently 
become apparent that the SFRs of normal star-forming galaxies 



are roughly proportional to galaxy stellar mass (e.g., Noeske et al.l 
| 2007l; lElbaz et al.|l200l iDaddi et alj|200l IPannella et alH 2009a: 
iDunne et alj|2009l ; Magdis et al.ll2010l) . meaning systematic differ- 
ences in the stellar masses of difference populations of galaxies 
(such as AGNs and non-AGNs) will also show as differences in 
their SFRs. Specific star-formation rates (i.e., the level of star for- 
mation per unit stellar mass; SSFR), by contrast, take this "mass ef- 
fect" into account, providing a measure of the relative growth rates 
of galaxies. As such, SSFRs give us a clear picture of how rapidly 
a given galaxy is forming stars relative to its star-formation history. 

Specific SFRs could also help to distinguish the fuelling mech- 
anisms that drive both nuclear and star formation activity. It is 
now becoming increasingly evident that almost all star formation 
takes place in one o f two regimes: the majority "main-sequence" 
(Noeske et ai] |2007h and the less common compact "starbursts", 
with the latter having significantly higher (i.e., factors of > 3) 
SSFRs. Interestingly, it seems that the type of regime is closely 
linked to the physical process driving the star formation activity, 
with starbursts being induced by major mergers, while internal pro- 
cesses (e.g., dis k instabilities, turbulence) drive main-sequence star 
format i on (e.g., lElbaz et ai]|2007l; iTacconi et alj|2008l ; IDaddi et al.l 
l2010bl ; iGenzel et al.l I2OI0I ; lElbaz et alj|201ll) . Since both major 
mergers and internal processes can als o channel gas and dust into 
the central regions of a galaxy (e.g. , Mihos & Hernquist 1 19961 ; 
iKormendv & Kenn icutt 2004l; lHopki ns & Hernquistl2006l) . it is rel- 
evant to question how nuclear activity fits within this regime. Until 
recently, nuclear activity was often considered to be closely linked 
to major merger events - largely due to the high fractions of highly 
luminous, "quasar class" AGNs that appear t o be associated with 
merg i ng systems (e.g., [Sanders et al.l Il988l; ISanders & Mirabell 
| 1996| ; ISuraceetal .11 19981 ; ICanalizo & Stocktonll200ll ; llvison et al.l 
120101) . However, a number of recent morphological studies of AGN 
hosts have presented evidence to suggest that more moderate levels 
of nuclear activity (i.e., sub-quasar; Lx^S 10 44 er g s s ~') are typi- 
cally associated with internally-evolving galaxi es undergoing sec- 
ular evolution rather than major-mergers (e.g ., [Grogin etalj |2005; 



ISchawinski et alj|201 ll ; ICisternas et alj|201 lh . On the other hand, 
measuring morphologies and identifying evidence of mergers can 
be problematic, especially at high redshifts where source faintness 
can introduce systematic effects. Specific SFRs, therefore, provide 
a direct, independent means to distinguish between merger and in- 
ternally induced star formation with which to verify such claims. 

In this paper, we use data from the deepest surveys yet un- 
dertaken at FIR wavelengths (i.e., GOODS-Herschel; hereafter, 
GOODS-H; P.I.: D. Elbaz) to measure the SFRs and SSFRs of X- 
ray selected AGN-hosting galaxies to z ~ 3. We use X-ray selec- 
tion as it provides the most accurate measure of nuclear activity, 
largely free from host contamination and less affected by obscura- 
tion than, for example, optical wavelengths. Throughout our anal- 
yses we use the Swift-BAT sample of AGNs as a local comparison 
sample as they have similar X-ray properties as the high redshift 
AGNs in the deep GOODS fields. Here, we largely focus on mod- 
erate luminosity AGNs (i.e., Lx= 10 42 -10 44 ergs s~') which are 
responsible for the majority of the integrated X-ray output of AGNs 
in the z < 3 Universe (e.g.JUeda et ai1l2003l ; lHasinger et al]|2005l; 
IHopkins et al.l2007l : lAird et al.l2010h . Although, whether they rep- 
resent the majority of black hole growth depends on the adopted 
bolometric correction factor^ Our principal aim is to compare the 



1 Here and throughout we use the term "nuclear activity" to explicitly refer 

to accretion onto a black hole (i.e., AGN) rather than, for example, intense 2 The luminosity-dependent correction factor of [Marconi et alj feOO- 



nuclear star formation. suggests that more luminous quasars dominate black hole growth, al- 
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SFRs and SSFRs of these moderate AGN hosts against the gen- 
eral galaxy population to establish whether (a) the levels of star 
formation activity around moderate luminosity AGNs differ from 
the general galaxy population and (b) merger or internally driven 
processes have dominated their stellar and black hole mass build 
up. 

This paper is structured as follows: in jj2]we describe our data 
acquisition and processing, while in Sj3]we describe the stacking 
procedure that we use to ensure we consider the FIR output of all 
X-ray AGNs in our sample, in addition to those that are formally 
detected. We describe how we convert FIR fluxes into SFRs and 
SSFRs in S|4]before presenting our results in S|5] where we demon- 
strate that moderate luminosity AGNs typically reside in main- 
sequence galaxies. In [j6]we interpret these results in terms of our 
current understanding of galaxy growth and consider the fraction of 
main-sequence galaxies containing AGNs. Finally, we summarise 
our findings in [jT] 



2 DATA 

Our main focus is the analysis of the FIR emission (measured by 
Herschel) of X-ray detected AGNs, as this provides a more ac- 
curate measure of t he SFRs of AGN hos ts than the mid-infrared 
(hereafter, MIR; e.g, Mull anev et alfeoiOh . However, an important 
aspect of our analyses also incorporates MIR data (from Spitzer) 
to derive the average MIR to FIR SEDs of AGNs in the Chan- 
dra Deep Fields (hereafter, CDFs; see i]5.3l >. We use optical and 
near-infrared (hereafter, NIR) photometry to estimate host galaxy 
stellar masses (hereafter, Af s tars)> which we combine with SFRs to 
calculate SSFRs. To aid in the interpretation of the deep (i.e., high 
redshift) CDF data, we make use of archival IRAS data for the Swift- 
BAT sample of local X-ray AGNs. In this section, we describe our 
X-ray AGN selection ( flOj the FIR, MIR, optical and NIR data 
and the cross-matching between these regimes ( i]2.2| and S\2.3\ . We 
end this section with a description of the Swift-BAT sample (i ]2.4l ). 



2.1 X-ray data 

Our initial sample selection is the 1243 X-ray sources detected 
in the 2 Ms CDF-N an d 4 Ms CDF-S surve ys, as reported in the 
lAlexander etafl J2003I) and lXue et all boill) main catalogues, re- 
spectively0 The vast majority of sources in these combined cata- 
logues have been identified as AGNs (i.e., 978, or ~79 per cent), 
with the remainder identified as either stars (i.e., 24, ~2 per cent) 
or starburst ing galaxies (i.e., 2 41, ~19 per cent). We use the classi- 
fications of Bau er et alj fc004h and lXue et al.l fcoill) to select only 
AGNs from the CDF-N and CDF-S catalogues, respectively. Both 
classification schemes combine information derived from the opti- 
cal (spectra and/or photometry) and X-ray regimes to differentiate 
between the different classes of object. Redshifts and absorption- 
corrected 2-10 keV luminosities (i.e., Lx) for each X-ray AGN 
in the CDF-N c atalog ue were taken from the catalogue described 
(2004). For the CDF-S we used the values from 
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Bauer et al 



IXue et aljfcoill) . converting their 0.5-8 keV to 2- 10 keV luminosi- 
ties by assuming a constant conversion factor, Lx= 0.68 io.5-8keV> 



though Vasudevan & Fabian 1 2007) argue against such a simple luminosity- 
dependent bolometric correction factor. 

3 We wish to make it clear that, throughout, we use the "CDFs" notation to 
refer to both Chandra Deep Fields and use CDF-S to refer to the Chandra 
Deep Field South only. 
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Figure 1. Intrinsic (i.e., absorption-corrected), rest-frame 2-10 keV lumi- 
nosities (i.e., Lx) of the X-ray AGNs detected by the 2Ms and 4Ms Chandra 
Deep Field Surveys (North and South, respectively) plotted as a function of 
redshift (see 32.11 . The AGNs that are detected (at > 3a) in either of the 
100 urn or 160 urn bands are circled (see key). To mitigate the effect of bi- 
ases when measuring the average infrared properties of X-ray AGNs in our 
sample we define a homogenous (in terms of Lx and z), restricted sample 
(shaded region). However, we stress caution when interpreting the results 
from the z = 2 — 3, Lx= 10 42 -10 43 ergs s bin (hatched region) due to the 
low number of X-ray detected AGNs in this redshift and Lx range (i.e., 29 
X-ray AGNs). 



appropriate for the constant spectral index (i.e., T = 1.8) used to 
derive their Lp s-8k aV values. I t shou ld be noted that the Lx es- 
timates reported by IXue et al. I Hh 1) were derived using simple 
band ratios to estimate the absorbing column densities, rather than 
the mo re detailed spectra l fitting appr oach carr i ed ou t by, for ex- 
ample, iBauer etHl, ( 2004). However, IXue et alj J20T 1) report that, 
when available, Lx values derived from spectral fitting generally 
agree (i.e., to within ~30 per cent) with those estimated using 
band ratios (although this factor may be higher for highly obscured 
sources). Of the 978 X-ray detected AGNs in the combined fields, 
879 (i.e., ~90 per cent) have either a spectroscopic (i.e., 480; ~55 
per cent) or photometric (i.e., 399; ~45 per cent) redshift measure- 
ment. The adopted X-ray luminosities have been corrected for ab- 
sorpti on on the basis of X-ray spectral analyses and X-ray band 
ratios dXue et al.ll2(5T 1; Bauer et al. in prep). As a key goal of this 
study is to explore the infrared properties of AGNs as a function 
of AGN power (i.e., Lx) and redshift, we only consider those 879 
AGNs for which these parameters are known. The Lx-z distribu- 
tion of these AGNs is shown in fig.Q] As for all flux limited surveys 
there is a clear Malmquist bias within this sample. Despite this, the 
regions of the parameter space spanning z = 0.5 — 3 and Lx= 10 42 - 
10 44 ergs s _1 are well sampled and we mainly focus our attention 
on the 564 AGNs within these ranges, although we stress caution 
when interpreting the results from the single z = 2 — 3, Lx= 10 42 - 
10 43 ergs s _1 bin, due to the low number of detected X-ray AGNs 
in this redshift and Lx range (i.e., 29 X-ray AGNs). 
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2.2 Mid and far-infrared data 

To characterise the MIR and FIR properties of our sample of X-ray 
AGNs, we used the deepest data currently available at 16, 24, 100 
and 1 60 \im: the Spitzer IRS peak-up survey (16 pm; Tepli tz et al.l 
1201 ll ; hereafter, GOODS-/i?5), the GOODS-Spitzer Legacy survey 
(24 urn; P.I.: Mark Dickinson), and the GOODS-H survey (100 tim 
and 160 mru lElbaz et alj|2~01 lh . The 16 yim observat ions and data 
reduction are described in full in lTeplitz et alj ( feoill) but, to sum- 
marise, they resulted in a pair of mosaics with pixel scales of 0.9" 
and average 3a depths of 40 and 65 pJy in the northern and south- 
ern fields, respectively. The 24 \im observations resulted in a pair 
of final mosaics with a pixel scale of 1.2" and a point-source 5a 
sensitivity limit of ~ 30 yJy in bot h fields. The GOODS-H observa- 
tions are described in Elba zetalJi20T J) and t he final images wer e 
produced following the procedure outlined in lBerta et alj f2oTob . 
In total, the GOODS-H program co nsists of 124.6 hours of PACS 
observations lPogl itschetafll2010l) in the northern field, and 264 
hours in the southern field (31.1 hours of SPIRE observations were 
also obtained for the northern field, but we do not use these longer 
wavelength data here). The final images have pixel scales of 1.2" 
(100 vim) and 2.4" (160 urn). At 100 pm the 5a sensitivities are 
1.7 mjy and 1.2 mjy in the northern and southern fields, respec- 
tively, while at 160 pm they are 4.5 mJy and 3.5-4.5 mjy, respec- 
tively0 

Although the various X-ray and infrared surveys are roughly 
centred on the same positions on the sky (i.e., [RA, Dec]: [12 :37 m , 
62°:14'] for the northern field; [3 /l :32 m , -27°:48'] for the southern 
field), because of the different fields of view and observing modes 
of the various instruments they do not cover exactly the same re- 
gion of the sky. Because of this, roughly 53, 70 and 47 per cent (i.e., 
465, 614, 409) of the X-ray AGNs in our sample with known red- 
shifts and L x are covered by the GOODS-IRS, GOODS-Spitzer and 
GOODS-H surveys, respectively. Roughly 44 per cent (i.e., 391) of 
our X-ray AGN sample is covered by all three surveys, although 
we only use data from all three when deriving the average infrared 
SEDs of the X-ray AGNs in our sample (see §5.31 ; we use only FIR 
fluxes from the GOODS-H survey to calculate Lir (see §4. It . 

To identify 16 pm counterparts to the X-ray AGNs in the 
CDFs we used the publicly available catalogue of Tepli tz et alj 
J201 1). The fluxes reported in the 16 pm catalogue were obtained 
by performing PSF-fitting using the positions of 5a sources from 
the 3.6 pm GOODS catalogue as priors (see Dickinson et al. in 
prep.). The 16 pm catalogue lists the fluxes of all sources detected 
at > 5a in the 16 pm band and contains 1309 sources across both 
CDFs. The 24 pm catalogue, which is also based on the positions 
of 5a, 3.6 pm priors, contains 5188 sources across both CDFs at 
> 3a. The 100 pm and 160 pm fluxes were calculated in a simi- 
lar fashion by performing PSF fitting at the positions of the > 3a 
24 pm priors. As we are exploring general trends in the FIR prop- 
erties of X-ray AGNs, rather than focusing on individual sources, 
we use a relatively low detection threshold of 3a at 100 pm and 
160 pm to maximise the number of sources in our sample. Across 
both GOODS-H fields there are 1324 and 1095 Herschel sources 
that satisfy this criteria at 100 pm and 160 pm, respectively. 

We used positional matching to identify infrared counterparts 
to the X-ray sources in each infrared band. Because there is no 
significant difference between the spatial resolutions of the 16 urn 
and 24 pm images we used the same 3.0" matching ra dius for both 
wavelengths, which we noted in lMullanev et al.l J2010) provides the 

4 The 100 pm and 160 pm surveys cover the same region of sky. 



best compromise between finding real matches while reducing the 
number of spurious matches. The 100 pm and 160 pm positions are 
from 24 pm priors (which are, in turn, based on IRAC 3.6 pm pri- 
ors), so there was no need to perform any matching for these bands. 
At 16 pm, we identified 192 matches (i.e., > 5a at 16 pm) to the 
X-ray detected AGNs with measured redshifts (i.e., ~41 per cent of 
those covered by this survey). At 24 pm we identify 449 matches 
(i.e., ~73 per cent) while at 100 pm and 160 pm we identify 174 
(i.e., ~43 per cent) and 140 (i.e., ~34 per cent) matches, respec- 
tively (i.e., > 3a in each of the 24, 100 and 160 pm bands). 0The 
larger numbers of matches at 24 pm compared to at 16 pm is due 
to the increased depth of the 24 pm survey. Approximately 25 per 
cent (i.e., 99) of the X-ray AGNs with measured redshifts covered 
by all three surveys are detected in all four bands. In each band the 
expected fraction of spurious matches is less than 2 per cent (cal- 
culate d using the same procedures as outlined in Mullan ey et al.1 
l20ld) . 

For each of the four infrared bands considered, the depth of the 
observations are different in the two fields. Furthermore, the depth 
of the observations also vary as a function of position within each 
field, especially in the 16 vim band. Therefore, when combining the 
individual fluxes to derive mean values, we take into account these 
differences by weighting as 1/error 2 of the measured flux. 

2.3 Optical and near-infrared data 

We use optical and NIR photometry to estimate the host galaxy 
masses of the X-ray AGNs in our sample (see §4.21 . A detailed 
description of the multi-wavelength catalog assembly and of the 
photometric redshift and stellar mass estimation will be given in a 
separate paper (Pannella et al. in prep.). In the following, and for 
the sake of completeness, we will briefly describe the data set and 
the procedures used. 

In the GOODS-N field we have built a PSF-matched multi- 
wavelength catalog with 10 passbands from U to 4.5 yim. The op- 
tical, NIR and IRAC data used here are presented in Cap ak et all 
(2004), W ang et all [2010) and Dickinson et al. (in prep.), respec- 
tively. A Ks-b and selected catalog has been built using SOURCE 
EXTRACTOR dBertin & Arnoutsl [l996T) in dual image mode. To 
compensate for the different resolution of the images we have ap- 
plied corrections based on the growth curve of "bona-fide" point- 
like sources. 

For the northern field photometric redshift determination 
clo sely follows the proced ure described in Pan nella et alj (2009b) 
an d lStrazzulloetalJfeOld). When comparing to the spectroscopic 
sample of iBarger et alj ( 120081) , we reached a relative (i.e., Az = 
kphot - Zspec|/[1+Zspec]) photo-z accuracy of 4 per cent with 3 per 
cent of catastrophic outliers (i.e. objects with Az > 0.2). For the 
southern field we have used the s pectroscopic or photometric red- 
shifts reported in lXue et alj J201lh thr oughout and u sed photome- 
try from the GOODS-MUSIC survey iSantini et ai]|2009l) to esti- 
mate ste llar masses (to ensure consistency, we fixed the redshifts to 
those o f IXue et alj|201 l] in these mass calculations). The GOODS- 
MUSIC catalog provides total magnitudes for a very similar pho- 
tometric cov e rage a s the GOODS-N field. We refer the reader to 
ISantini et alj J2009h for a detailed description of the catalog con- 
tent. 

5 Here, the percentages refer to the fraction of sources detected in the re- 
gion covered by each individual survey, i.e., not the fraction covered by all 
three surveys. 
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We used positional matching to identify the counterparts to the 
X-ray AGNs in the above catalogues. Rather than matching directly 
to the X-ray positions themselves (which have a high positional 
uncertainty, especially at large off-axis angles ), we used the posi- 
tion s of their optical counterparts reported in iBarger et alj J2003I) 
and lXue et al.l d201lh for the northern and southern fields, respec- 
tively. We use a matching radius of 0.5"in both fields and match 
?«95 per cent (i.e., 388) of the 409 X-ray AGNs with redshift mea- 
surements and GOODS-H coverage. 

We found that for 27 (i.e., ~ 7 per cent of the total sample cov- 
ere d by GOODS-H) A GNs in the northern field the redshift from 
the IBarger et al.1 J2003h differed from our own new estimate by a 
significant amount (i.e., Az > 0.2). We r eject these sources en tirely 
from our analyses as we believe that the IBarger et al.l d2003l) pho- 
tometric redshift is now out-of-date (i.e., has been superseded due 
to improved photometric/spectroscopic data) meaning the value of 
Lx is likely incorrect. This ensures that, across both fields, we have 
used fully consistent redshifts throughout to derive Lx, Lir and 
M stars . 



2.4 The Swift-BAT sample of z ~ AGNs 

As we described in iMullanev et alj i l2010h . the Swift-BAT sample 
of AGNs jTueller et all l2008h iscurrenflv one of the best local 
comparison samples for studies of high redshift, X-ray AGN s de- 
tected in the CDFs (see Figs. 7 and 10 of Iwinter et ai]2009l where 
they plot the L x and iV H of the Swift-BAT AGN sample). Of the 
104 Swift-BAT A GNs for which Lx and Nu have been measured 
dWinter et al]|2009l) . 41 have infrared flux density measurements 
from all four IRAS bands (these 41 objects are hereafter referred 
to as the BAT/IRAS sample), allowing direct comparison with the 
infrared properties of the AGNs in the CDFs. All of these AGNs 
have z < 0. 1 and are thus "local" in comparison to the z = 0.5-3 
CDF AGNs which form ou r high redshift sample. As none of the 
Bl azars in the Tueller et alj q2008h have Lx measurement reported 
in Iwinter et al.1 J2009h they are automatically excluded from our 
sample. 



3 MID AND FAR-INFRARED STACKING PROCEDURE 

By using the deepest FIR observations currently available, we de- 
tect roughly ~46 per cent of X-ray sources at > 3a in either 
the 100 \im or 160 |i.m bands. This represents a vast improve- 
ment over previous F IR surveys (e.g., ~ 1 1 and ~ 20 per cent for 
Mull anev et aT]|2010l and lShao et alB oiO). respectively. However, 
we resort to stacking techniques (see i|5.1.2| > to confirm that the 
trends observed in the FIR-detected fraction of our sample remain 
when all X-ray AGNs in our sample are considered (i.e., even those 
undetected at FIR wavelengths). 

As our intention is to explore how the infrared properties (e.g., 
infrared luminosities, SED shape) change as a function of redshift 
and nuclear activity, we stack subsamples of our X-ray AGNs sepa- 
rated into six different Lx and z bins: all combinations of z =0.5-1, 
1-2, 2-3 and L x = 10 42 -10 43 ergs s^'and 10 43 -10 44 ergs s _1 ; the 
number of AGNs in each stack are given in column 4 of table [T] 
These bins cover the region of the Lx~z parameter space that bene- 
fits from high levels of completeness (see ^2. 1 l and fig.[U(. The na- 
ture of small, yet deep fields such as the CDFs implies that, while 
we can probe moderate luminosity (i.e., Lx= 10 42 -10 44 ergs s _1 ) 
AGNs to high redshifts, there are only a small number of the rarer, 



quasar luminosity AGNs in our sample, especially at lower red- 
shifts (i.e., z S, 2). W ide, but less deep, infrared surveys such as 
PEP jLutz et alJ[20lTh and HERMES (PL: S. Oliver) are needed 
to fully explore the average infrared properties of both higher lu- 
minosity (i.e., Lx> 10 44 ergs s ) or, alternatively, lower redshift 
(i.e., z < 0.5) AGNs. 

We only stack on residual images (i.e., images from which the 
sources down to 3a have been removed), stacking only those X-ray 
sources that are undetected in the various infrared bands (i.e., < 5a 
for 16 vim and < 3a for 24 |im, 100 \im and 160 |i.m). We take 
this approach, rather than stacking all X-ray positions on the origi- 
nal images, to reduce the effect of flux boosting from nearby bright 
sources. However, as the residuals from very bright sources can be 
strong even after source subtraction, we do not stack at positions 
that lie within 10" of > 5a sources (equivalently, >65— 1 10 (oJy, 
> 1 .2—1 .7 mjy and >4-4.5 mjy at 16 [im, 100 \im and 160 \im, re- 
spectively). The stacked images are produced by summing cut-out 
images centred on the positions of the infrared-undetected X-ray 
sources in both fields. Prior to summation, each pixel in the indi- 
vidual cutouts is weighted by the inverse of its RMS noise, deter- 
mined from RMS maps of the fields. The weighted mean stack is 
then calculated by dividing the sum of these cutouts by the sum of 
the weights on a pixel-by-pixel basis, i.e., 

L?=iSi(x,y)Wi(x,y) 



S(x,y) = 



rUMx,y) 



(1) 



where S(X,y) is the final stacked image, 5,(x,v) and w,-(x,y)(= 
l/RMS 2 [xj]) represent the individual cutout images and weight 
map, respectively, and n is the total number of cutouts to be com- 
bined in the stack. Once the weighted mean stacks have been pro- 
duced for the sources in each of our bins, we use aperture photom- 
etry to measure the flux within a circular aperture placed at the cen- 
tre of each stacked image. The radius of this aperture depends on 
the PSF of the observation and thus, ultimately, on the wavelength 
and instrument: 6" for both the 16 |J.m and 24 \xm wavebands and 
9" and 15" at 100 |J.m and 160 (im, respectively. Through simu- 
lated stacks, we find that the 100 and 160 urn fluxes extracted 
through aperture photometry are 70 per cent of those derived from 
PSF fitting (i.e., the technique used to extract individually detected 
sources). To correct for this discrepancy, we apply an aperture cor- 
rection to our stacked fluxes, such that the reported flux is 0.7 -1 
times the flux derived from the stacks. We also include a 15 per cent 
upwards correction factor to account for the effects of source mask- 
ing that was performed when removing the low-frequency noise 
from the raw Herschel images. We then use a weighted mean (here, 
weighting by number) to combine this flux derived from stacking 
with the mean flux from the individually detected sources. We fol- 
low the same procedure (aside from minor modifications to account 
for the different image resolutions) to stack in each of the four in- 
frared bands. 

The errors on the fluxes derived from these stacks of the unde- 
tected sources are calculated by using the same procedure to stack 
the same number of random positions in the residual frames. We 
perform this "random stacking" 10,000 times, with the error cal- 
culated by taking the standard deviation of the fluxes from these 
trials. 

We use a different bootstrapping approach to calculate the er- 
rors on the total mean 100 \im and 160 p_m fluxes (i.e., those de- 
rived by combining the detected and undetected mean fluxes). This 
"error on the mean" is different from both the error on the stacked 
fluxes and the standard deviation of the samples. Taking each bin 
in turn, we randomly select (with replacement) 1 of the sample 
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Figure 2. Observed [100 u.m/24 urn] (left) and [160 u.m/24 p] (right) flux 
ratios vs. source redshift for the X-ray AGNs in our sample. Blue squares 
and red circles represent Lx< 10 43 ergs s _1 and > 10 43 ergs s AGNs, 
respectively. Downward pointing arrows indicate 3cj upper limits for those 
X-ray AGNs that are undetected at FIR wavelengths. We use the observed 
100 [im and 160 |xm fluxes to measure Ltr for AGNs at z = 0.5 — 1.5 and 
Z = 1.5 — 3, respectively (see text, §4. Q . Included in these plots are tracks 
indicating the expected flux ratios if 25, 50 and 100 per cent of the flux in 
the 100 [im and 160 iim b ands arises from the A GN (estimated using the 
empirical SED templates of Mullaney et al. 201 1). We exclude from further 
consideration any systems in which these ratios indicate that more than 50 
per cent of the 100 |a.m or 160 flux (depending on redshift) is dominated 
by the AGN (16 in total, see text; igj). 



within that bin, taking no account of whether they are detected or 
not, and use the procedure outlined above to calculate the weighted 
mean flux of this sample. This is performed 10,000 times for 
each bin and the error on the mean is calculated by taking yVJ 
times the standard deviation of the results of these trials. 



4 CALCULATING TOTAL INFRARED LUMINOSITIES, 
STAR FORMATION RATES, HOST GALAXY MASSES 
AND SPECIFIC STAR FORMATION RATES 

One of the main aims of this study is to relate nuclear activity to star 
formation at z < 3. A key measurement is therefore the SFRs and 
SSFRs of the AGN hosts in our sample. In this section we describe 
how we measure Lir for the AGN hosts in our sample and convert 
these into star formation rates. We also describe how we use deep 
optical and NIR observations to measure the stellar masses of these 
galaxies, enabling us to calculate their SSFRs. 



4.1 Total infrared luminosities and star formation rates 



We use the SED library of ICharv & Elbaj feOOll) to calculate the 
integrated 8-1000 [im infrared luminosities (hereafter, Lir) of the 
AGN hosts in our sample. These SEDs, when combined with the 
monochromatic luminosities measured by the 100 [im and 160 [im 
Herschel-VACS bands, have been shown to provide reliable es- 
timat es of Lir at the reds hifts considere d here (e.g., lElbaz et all 
20101) . The shapes of the ICharv & Elba j fcOOll) SEDs are purely 
a function of infrared luminosity, meaning that each SED has a 
single characteristic Lir (and vice-versa). We identify the SED 
that most closely matches the rest-frame monochromatic luminos- 
ity (i.e., vLy) derived from either the observed 100 [im or 160 (xm 
flux then integrate this SED between 8-1000 (xm. The choice of 
FIR waveband depends on the redshift of the source; for z =0.5-1.5 
sources, we use the 100 urn band, while for z = 1 .5-3.0 sources, we 
use the 160 [im band. This ensures that we always sample roughly 
the same rest-frame wavelengths (i.e., 40-67 [im), thus minimising 



the systematic uncertainties introduced by -corrections. Roughly 
~42 per cent (i.e., 144) of the z = 0.5 — 3 X-ray AGNs in our sam- 
ple that are covered by the GOODS-H footprint are detected in their 
respective FIR bands; hereafter, we collectively refer to these as be- 
ing "FIR detected". 

As the Charv & Elbaz (2001) catalogue only includes tem- 
plates for non-AGN-hosting galaxies, we must ensure that these 
rest-frame wavelengths are not dominated by the AGN, otherwise 
the Lir of the host galaxies will be overestimated. For this, we cal- 
culate the observed-frame [100 (im/24 p.m] and [160 iim/24 \im] 
flux ratios for the z =0.5-1.5 and 1.5-3.0 AGNs, respectively, de- 
tected in these bands. These ratios are presented in the left and right 
panels of fig. [2] Included in these plots are tracks indicating the ex- 
pected ratios if 25, 50 and 100 per cent of the flux in the 100 \im 
and 160 vim bands arises from the AGN. These tracks are calcu- 
lated by combini ng the empiri c al ave rage-AGN and "SB5" SED 
templates from lMullanev et alj ( 1201 II) and are appropriate for the 
Lx and Lir ranges considered here. We only calculate Lir when 
at least 50 per cent of the 100 p.m or 160 [im flux (depending on 
redshift) is attributable to the host galaxy. The tracks shown in fig. 
[2] are largely unaffected by emission or absorption features in the 
host galaxy templates because the 24 p.m band is always dominated 
(i.e., > 70 per cent) by the AGN when < 75 per cent of the observed 
100 |i.m and 160 |i.m flux arises from the host galaxy. That said, the 
levels of extinction applied to the AGN c omponent can a l ter the se 
tracks. When the average AGN SED from Mull aney et al] J20T ID is 
used (which already includes a modest amount of extinction) we re- 
move a total of 16 FIR-detected X-ray AGNs from our high redshift 
(i.e., CDF) sample (i.e., ~6 per cent of the sample). If, instead, we 
were to use a heavily extinguishe d AGN component (su ch as that of 
ESO103-035 orNGC 5506; see lMullanev et alj|201 ]]) the number 
of rejected sourced would increase by 18 (i.e., by a factor of ~ 2), 
but this would not affect our overall conclusions. We identify 16 of 
the 41 BAT/IRAS AGNs (i.e., ?»40 per cent) with 60 \wi flux likely 
dominated by the AGN, rather than the host galaxy. Since such a 
high fraction of the BAT/IRAS comparison sample is found to be 
AGN dominated, removing them entirely from our analyses would 
give a false impression of the average Lir and star-formation ac- 
tivity of local AGNs (especially at Lx > 1 43 ergs s _ 1 , where two- 
thirds of BAT/IRAS AGNs are "AGN dominated"). Instead, we as- 
sume that their rest-frame 60 \im fluxes set a hard upper limit on the 
infrared luminosities of their host galaxies. Star formation rates are 
calculated from the host galaxy infrare d luminositi es using the pre- 
scription outlined in Kennicutt (1998) assuming a lSalpeteJ fl955T) 
IMF, i.e., 



^5 rr = 1.7xl0- 10 (^ 



M Q yr 



(2) 



3.8 x 10 33 ergs s _1 . 



and assuming L Q = 

From the rest-frame ultraviolet luminosities we estimate that 
incl uding the contribu tion from unobscured star formation (as in 
e.g jDaddi et alj2007l) would have a negligible effect on our results 
(i.e., at the level of l.Oln 7 P er cent) and, for simplicity, is neglected 
from the rest of our analyses. 

We use the same procedure to calculate the mean Lir and 
SFRs of the AGN hosts by treating the mean 100 pin or 160 (im 
flux (derived from stacking analyses) as an individual source at the 
mean redshift of the subsample. Errors on these average Lir and 
SFRs are calculated by passing the average FIR flux ± lo (see ^3} 
through this process. 



GOODS-H: Star formation in AGN host galaxies 1 




Figure 3. Left: Rest-frame U — V vs. V — J colours for the AGNs in our sample, over-plotted onto the colour-colour distribution of non-AGN galaxies in 
the same fields (grey points). We use contours to show the density of the non-AGNs in this plot. In general, type 2 AGNs tend to be redder than the galaxy 
population, which we argue indicates that their rest-frame optical to NIR SEDs (which we use to estimate host galaxy stellar masses, M slars ) are dominated by 
the host galaxy, rather than the AGN (see ^4.21 . As a precaution, we only use host galaxy stellar masses in cases where U — V > 1.0 and V — J > 1.0 (indicated 
by the dashed lines). Right: Histogram showing the distribution of V — J colours for both Type 1 and Type 2 AGNs in our sample and the underlying galaxy 
population. For clarity we have expanded vertical scale of the galaxy fraction by a factor of 3. 
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Figure 4. Example Hubble Space Telescope ACS three-colour (i.e., B, V, I) images of the galaxies identified as having "secure" stellar mass measurements 
(see EI4.21 . We primarily use colour-colour selection to identify those galaxies whose integrated light is likely dominated by the host-galaxy (and therefore 
gives a reliable measure of the stellar mass of the host galaxy). As confirmation of this, we note that > 90 per cent of z < 1 galaxies selected in this way also 
show clear evidence of having extended morphologies and/or are not dominated by a central, point-like source (i.e., AGN). At higher redshifts it becomes 
increasingly difficult to assess galaxy morphologies using optical imaging, so we rely solely on rest-frame colour-colour selection at z > 1 . The diameter of 
the small circle in the top-right of each image corresponds to the sp atial resolution o f the i nstrument i n the /-band. The labels are a combination of the field 
(N=CDFN, S=CDFS) and the index or the X-ray source from either lAlexander et alj fe003l: CDFN) or lXue et alj feOllI : CDFS). The crosshairs indicate the 
position of the optical counterpart to the X-ray source. 



4.2 Host galaxy stellar masses and specific star formation 
rates 

Stellar masses were estimated for both fields us ing t he SEP fit- 
ting co de described in detail in lDrorv et al.l d2004l) and lDrorv et al.l 
( 2009). We parametrize the possible star formation histories (SFHs) 
by a two-component model, consisting of a main, smooth compo- 
nent described by an exponentially declining star formation rate 
\|/(f) °c exp(— f/t), linearly combined with a secondary burst of 
star formation. The main component timescale x varies in £ 
[0.1,20] Gyr, and its metallicity is fixed to solar. The age of the 
main component, t, is allowed to vary between 0.01 Gyr and the 
age of the Universe at the object's redshift. The secondary burst of 



star formation, which cannot contain more than 10% of the galaxy's 
total stellar mass, is modelled as a 100 Myr old c onstant s t ar for - 
mation rate episode of solar metallicity. We adopt a ISalpetenfl955f) 
initial mass function for both components, with lower and upper 
mass cutoffs of 0.1 and 100 M Q . 

Adopting the ICalzetti et al] (2000) extinction law, both the 
main component and the burst are allowed to exhibit a variable 
amount of attenuation by dust with Ay 2 £ [0, 1 .5] and [0, 2] for 
the main component and the burst, respectively. For the GOOPS-N 
field, we have corrected the "aperture magnitude" stellar masses to 
total stellar masses by means of the ratio of total (FLUX_AUTO) 
and aperture fluxes in the A"-band detection image. Of the 346 host- 
dominated X-ray AGNs with GOOPS-H coverage and consistent 
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redshifts (see t]2.3t . we have been able to calculate stellar masses 
for 327 (i.e., ~95 per cent). 

Of course, our mass estimates are based on the assumption that 
the light in these optical and NIR bands is dominated by the host 
galaxy, rather than the AGN. While a number of previous stud- 
ies have found that this is indeed the case forL x < 10 44 ergs s 
AGNs such as those considered h ere (e.g. Isilverman et al.lr2008l : 
iBundv et alj2008l : IXue et alj20ld) . we wanted to perform our own 
tests to confirm that this is true for our sample. 

To test whether the optical/NIR fluxes of our AGNs are host- 
dominated we compare their rest-frame colours against those of 
the general galaxy population. The central engines of AGNs have 
blue optical-NIR colours as their broad-band UV to NIR SEDs rise 
monotonically to shorter wavelength - a result of this light being 
emitted directly from the accretion disk itself. Therefore, the pres- 
ence of a dominating AGN will tend to result in the galaxy as a 
whole (i.e., stellar + AGN components) having a bluer colour than 
if the AGN were not present. In the left-hand panel of fig.[3]we plot 
the rest-frame U — V and V — J colours of the AGNs in our sam- 
ple, separated into Type 1 (i.e., optically unobscured; broad-line) 
and Type 2 (i.e., optically obscured; narrow-line) AGNs0 Also 
included in this plot are the rest-frame colours of the underlying 
galaxy population in the GOODS fields. Rather than having par- 
ticularly blue colours, Type 2 AGNs tend to lie among the redder 
galaxies in this colour-colour plot, indicating that their active nu- 
clei are obscured and that their optical-NIR SEDs are likely host- 
dominated. Furthermore, there is no evidence to suggest that more 
X-ray luminous Type 2 AGNs have systematically different U — V 
and V — J colours. This is more clearly shown in the right-hand 
panel of fig. [3] where we plot histograms of the V — J colours 
of Type 1 and Type 2 (separated according to Lx) AGNs and the 
underlying galaxy population. Type 1 galaxies, on the other hand, 
tend to have much bluer colours, as may be expected if the unob- 
scured light from the accretion disk is dominating the output of the 
total galaxy. As a precaution, we only use galaxy stellar masses 
for those AGNs in our sample with rest-frame U — V > 1.0 and 
V -J > 1.0. This leaves 255 AGNs in the GOODS-H fields with 
"secure" masses, of which 1 15 are FIR detected. However, we note 
that our conclusions regarding the host galaxy stellar masses and 
SSFRs do not change if we, instead, include all AGNs in our anal- 
yses. We note that > 94 per cent of the AGNs in our sample that 
satisfy these criteria also satisfy \li i L] i NxLx> 2, which lBluck et al.l 
d201ll) suggest is a reliable indicator of a host-dominated optical- 
NIR SED. As a final test, we exploit Hubble Space Telescope B, 
V, I, z (i.e., ACS) images to check the optical morphologies of the 
selected galaxies to ensure that they are not dominated by a bright, 
central point source (i.e., AGN; see Ammons et al. 2 01 ll for a de- 
tailed discussion on this point) Q We find that > 90 per cent of z < 1 
galaxies in our sample have clear, extended morphologies (see fig. 
[4^, indicating that their integrated light is dominated by the host 
galaxy, rather than the AGN. At higher redshifts, it becomes in- 
creasingly difficult to assess optical morphologies, so we assume 
that our rest-frame colour cuts provide the same success rate at 
all redshifts. However, we also note that recent studies exploiting 
new NIR data from WFCAM3 on-board Hubble (i.e., the CAN- 
DELS survey; RI.s: H. Ferguson, S. Faber) find that AGN-hosts 



Rest-frame colours are calculated in the lBuser & Kurucj jl978t) U and V 
niters and in the / band transmission curve of the WIRCam instrument. 
7 Images were retrieved via the ACS GOODS Cutout Service v2.0: 
http://archive.stsci.edu/eidol_v2.php 



at these high redshifts have quantitatively similar morphologies to 
their non-AGN counterparts, again suggesting that their integrated 
light is dominated by the host galaxy (D. Kocevski, private com- 
munication). 

We calculate mean stellar masses for each of our Lx and z 
bins and use bootstrapping to calculate the error on these means in 
a manner consistent with that use to determine the errors on our 
mean 100 ixm and 160 ii.m fluxes (and hence, Lir and SFRs; see 
ij3}. Taking each bin in turn, we randomly select (with replacement) 
/3 of the sample within that bin, and calculate the mean stellar 
mass of this 1 sample. This is performed 10,000 times for each 
bin and the error on the mean stellar mass is calculated by taking 
\/ 1 times the standard deviation of the results of these trials. 
Errors on the SSFRs are derived by combining the errors on the 
SFRs and stellar masses in quadrature. Mean stellar masses and 
their associated errors for the AGNs in each of our Lx and redshift 
bins are given in tableQ] 



5 RESULTS 

In this section we present our results and explore the links between 
nuclear activity, measured from Lx, and the levels of star formation 
activity, measured from Lir, taking place in their host galaxies. We 
also pay particular attention to how star formation activity in AGN 
host galaxies changes with redshift. First, we consider the ~ 40 per 
cent of X-ray AGNs for which we are able to calculate host galaxy 
Ljr values and SFRs, noting that we have already excluded any 
galaxies where the AGN likely dominates the observed FIR flux 
(see ^4. It . Then, we use the results from our stacking analyses to 
confirm that trends seen in this FIR-detected fraction of our sample 
extend to the general X-ray AGN population. 

5.1 Infrared luminosities and star formation rates of 
AGN-hosting galaxies 

5. 1. 1 X-ray AGNs detected at 100 /jm and 160 yam 

In the left-hand panel of fig.|5]we plot Lir versus Lx for all the FIR- 
detected AGNs in both the high redshift CDFs and local BAT//&4S 
samples. Over the entire sample, which spans over three orders of 
magnitude in both Lir and Lx, we find no statistically significant 
correlation between Lir and Lx (Spearman's rank correlation co- 
efficient: 0.29). Furthermore, this Lir-Lx correlation typically gets 
even weaker when we consider the CDF AGNs in discrete redshift 
bins (i.e., z =0.5-1, 1-2, 2-3; Spearman's rank: 0.16, 0.09 and 0.30, 
respectively). We note that there is a slightly stronger Lir-Lx cor- 
relation within the BAT/IRAS sample (Spearman's rank: 0.36) due 
to the rise in Lir of the BAT sample at Lx> 10 43 ergs s _1 . How- 
ever, this is likely due to the large fraction of the AGNs in this low 
redshift sample BAT/IRAS being AGN dominated at 60 \im (see 
S|4.1b . Indeed, the minimum Li r produced by a "pu re" AGN in- 
creases almost linearly with Lx dMullanev et al.l201ll : indicated by 
the dotted line in fig. [5}. By splitting the sample into redshift bins, 
we find that, generally, the Lir of AGNs increase strongly with 
redshift. Indeed, the AGNs in our highest redshift bin (i.e., z =2-3) 
typically have Lir and, correspondingly, SFRs around 1.5 orders of 
magnitude higher than the z < 0. 1 BAT/IRAS AGNs with the same 
X-ray luminosities (see ^5. 1 .21 for a more quantitative breakdown 
of the increase in the average Lir with redshift). The weak cor- 
relation in the unbinned sample arises from a combination of this 
strong Lir-z correlation and an increase in the proportions of more 
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Figure 5. Left: Total infrared luminosities of X-ray AGNs vs. intrinsic (i.e., absorption-corrected) rest-frame 2-10 keV luminosities (i.e., Lx) for the far-infrared 
detected X-ray AGNs in our sample. We have separated the sample in terms of redshift (see key). Open stars are used to indicate those AGNs in the local 
(i.e., z < 0.1) BAT/ffiAS comparison sample whose 60 u.m flux is likely dominated by emission from the AGN, rather than the host galaxy. In such cases the 
Lir shown here can be considered a hard upper limit on the infrared luminosity of the host galaxy (see ij4. 1 1 . In a given redshift bin there is no significant 
correlation between X-ray and infrared luminosities, but we see a strong increase in infrared luminosity with increasing redshift. Right: Same as left, but here 
we show the average infrared luminosities derived from stacking analyses (see text, [[5J vs. mean-average Lx for each of our Lx and redshift bins. Since 
two-thirds of L x = 10 43 -10 44 ergs s AGNs from the BAT/IRAS comparison sample only have upper limits on their host galaxy Lir, we estimate an upper 
limit for the average value. Again, we see a clear increase in Lir with redshift, but no significant change over the L x = 10 42 - 10 44 ergs s 1 range considered 
here. In both panels w e indicate the Lir expected from a "pure" AGN (i.e., without any host-galaxy contamination; dotted line), calculated using Eq. 4 of 
iMullanevetaklfcOllI) . 



X-ray luminous AGNs in our sample at higher redshifts; a result of 
both selection effects and a real , observed increase in the average 
Lx w ith increasing redshift (e.g., Barg er et alj2005l ; lHasinger et al.l 
l2005t) . 

The increase in Lir with increasing redshift is more clearly 
visible in fig. [6] where we plot Lir as a function of redshift. Here, 
we have included upper limits on Lir, calculated in the same way 
as for the detected sources assuming 3a upper limits on the FIR 
fluxes. Qualitatively, the distribution of AGNs in the Lir-z plane is 
the same as that of the FIR-detected non-AGNs in the same field. 
As is expected when ~60 per cent of a sample is not detected, this 
plot shows a clear Malmquist bias. However, it is clear that the 
upper envelope of the Lir-z relation increases strongly with red- 
shift; by about 1.5 orders of magnitude from z = 0.5 to 3.0, i.e., 
Lir~2x10 u to 10 13 Lq. This trend continues to the low redshifts 
(i.e., z < 0.1) spanned by the BAT sample and is too steep to be 
explained in terms of the increase in the comoving volume of the 
survey expanding the range of Lir at higher redshifts; the probed 
comoving volume between z = 2 — 3 is only a factor of 4 times 
larger than that between z = 0.5 — 1. 

We have demonstrated that the FIR-detected AGNs in both the 
high redshift CDF and low redshift BAT samples show no strong 
increase in Lir with Lx, yet are consistent with a strong increase 
in Lir with increasing redshift. We will look more closely into the 
causes and consequences of these results in i ]5.2l and fj6]but, first, 
we will confirm whether these results still hold when we consider 
all X-ray AGNs in the GOODS fields, i.e., not just those detected 
at FIR wavelengths. For this, we resort to stacking analyses. 



5.1.2 Results from stacking analyses 

By considering only the detected sources in the previous subsec- 
tion, we sampled only the most intrinsically infrared luminous 
galaxies at each redshift. In this subsection we consider the con- 
tribution of X-ray AGNs that are not formally detected in each of 
the wavebands by relying on stacking analyses (see ij3} to push the 
average detection threshold beyond what is achievable for individ- 
ual sources. 

In the right-hand panel of fig. [5] we plot the average Lir (de- 
rived from stacking analyses) vs. average Lx for CDF AGNs in 
each of our two Lx bins and three redshift bins (i.e., Lx = 10 42 - 
10 43 ergss" 1 and 10 42 -10 43 ergss" 1 ; z =0.5-1, 1-2 and 2-3). 
Also included in this plot is the mean-average Lir vs. mean-average 
Lx of the low redshift BAT /IRAS AGNs. These averages reproduce 
the same overall trends as those seen in the FIR detected fraction 
of our sample, i.e., that, within each redshift bin there is no strong 
correlation between the average Lir and average Lx, but that the 
average Lir increases strongly with increasing redshift. The aver- 
age properties, including the average Lir, of the AGNs in each of 
our Lx and redshift bins are presented in tableQ] 

The increase in the average Lir of our AGNs with increasing 
redshift is clearly seen in the right-hand panel of fig. [6] where we 
plot the average Lir as a function of redshift. The mean Lir of 
Lx=10 42 -10 43 ergs s -1 AGNs increases by a factor of 3. 5lJ'2 from 
z = 0.5 - 1 to z = 1 - 2 and by a factor of 10.3+^f, from z = 0.5 - 1 
to z — 2 — 3. The large uncertainty on the latter of these ratios is 
due to the small number of Lx=10 42 -10 43 ergs s AGNs in our 
highest redshift bin (i.e., 29, of which only 11 are detected at FIR 
wavelengths). Comparison with the Lx=10 42 -10 43 ergs s _1 AGNs 
in the z < 0.1 BAT/IRAS sample (mean Lir=2.23+o;^ x 10 10 L Q ) 
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Figure 6. Left: Rest-frame infrared luminosities vs. redshift for our sample of X-ray AGNs, split in terms of intrinsic, rest-frame 2-10 keV luminosities (i.e., 
Lx; see key). Downward pointing are used to indicate the 3a upper on Lir in cases where there is no formal FIR-detection. The points around z~0 represent 
the BAT//RA5 sample of local AGNs, which cover roughly the same range of X-ray luminosities as the high redshift AGNs in the Chandra deep fields. While 
there is a clear Malmquist bias, the increase in the upper envelope of the distribution with redshift is too large to be explained in terms of the increase in the 
surveyed comoving volume with increasing redshift (see j]5.1.U . The distribution of AGNs in this plane is similar to those of FIR-detected non-AGN hosts 
(shown as grey points in both panels) Right: The same as the left panel, but now showing the average Lir (derived from our stacking analyses; see J3) and Lx 
values. Again, the increase in Lir, and therefore host galaxy star formation rates (shown by the right hand axes), is clear (see i|5.1.2) . 



reveals an increase in the mean Lir by factors of 5.3 + | 

and 55±29 from z < 0.1 to z = 0.5- 1, z = 1 -2 and z = 2-3, 

respectively^ 

We see a similar increase in Lir with redshift for the more X- 
ray luminous AGNs in our sample (i.e., Lx= 10 43 — 10 44 ergs s ). 
Considering just the high redshift (i.e., CDF) sample first, there 
is a factor of 3.5 + [ '\ increase in the mean Lir between the z = 
0.5 — 1 and z = 1 — 2 bins and a factor 11.3^1'2 increase be- 
tween the z = 0.5 — 1 and z = 2 — 3 bins. Comparing with the 
mean Lir of BAT/IRAS sample AGNs from the same Lx interval 
(Lir=4.5+ 3 ; 3 ) x 10 10 L ) reveals an increase in Lir by factors of 
at least 2.9^% 10+" and 33 +3 j from z < 0.1 to z = 0.5 - 1, 
z = 1 — 2 and z = 2 — 3, respectively^ A similar increase in the 
Li r of L x = 1 43 -10 44 er gs s' 1 AGN s with redshift was reported 
by lLutz etal.1 ( l201Cj) and ishao et al] fcOld) . These observed in- 
cr eases in Lir are also consistent with the upper limits reported 
in lMullanev et al] J20T0T) that were derived from Spitzer 70 (i.m ob- 
servations. 

As a check, we also calculated the mean Lir of AGNs in the 
full Lx= 10 42 — 10 44 ergs s _1 bin and find that they are consistent 
with the results from the narrower Lx bins (see fig.[6]and table[T](. 
A least-squares fit to the average Lir of these AGNs gives: 

L :R /L = 1.4 x 10 13 (r/Gyr)- 24 (3) 

where t is the age of the Universe in Gyr. The influence of the FIR- 
undetected AGNs on these average Lir values is clearly seen in the 

8 The errors on the ratios presented in this section reflect the maximum 
and minimum possible values, the statistical errors will be roughly 30-40 
per cent smaller. 

9 In reality these ratios are likely considerably higher since the majority 
of L x = 10 43 -10 44 ergs s~ 1 AGNs from the BAT//RA5 comparison sample 
have only upper limits for Lir. 



z = 2 — 3 bin, which lies close to the bottom of the distribution of 
FIR-detected non-AGNs in the same field. This explains why the 
average Lir of our sample does not rise as quickly with redshift as 
the FIR-detected galaxies. 

Finally, we note that if we instead use volume-averaged bins, 
(=< Lir > A'agn/V(zi,Z2)> where Nagn is the number of AGNs 
in each redshift bin and V (21,22) IS the surveyed volume between 
the up per and lower redshift bounds, zi and Z2\ e.g jHeckman et all 
2004) we find that the infrared luminosity density of Lx= 10 42 — 
10 44 ergs s" 1 AGNs increases by a factor of ~3 from z=0.5-l to 
z=2-3. This is in general agreement with the incre ase in the to- 
tal SFR density over the s ame redshift range (e.g.. Ilvlurphv et ail 
2009 ; [Magnellietaitolll) and implies that the relative galaxy to 
black hole growth in moderate luminosity AGNs was higher at high 
redshifts compared to today (assuming that the bolometric correc- 
tions, accretion histories etc. have, on average, remained constant 
throughout; see also t\6.2\ . 

We have shown that the GOODS-H observations confirm the 
results from previous studies of the (far) infrared properties of X- 
ray selected AGNs over large redshift ranges; i.e., that, within the 
L x = 10 42 -10 44 ergs s _1 interval there is no significant L x -Lir cor- 
relation but Lir and thus SFRs, increases strongly with redshift. In 
the following subsection we explore whether this increase in Lir 
with redshift is caused by changes in the stellar mass or the SSFRs 
(or a mixture of both) of the AGN hosts. We find that it is most 
likely due to the latter and is broadly consistent with the increase 
in the SSFRs of star-forming galaxies with increasing redshift. 



5.2 Masses and specific star formation rates of AGN-hosting 
galaxies. 

In their study of the sub-mm properties of X-ray selected AGNs, 
iLutz et al] d2010t) showed that the increase in the sub-mm lumi- 
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Figure 7. Left: Host galaxy stellar masses (i.e., M stius ) vs. redshift for the FIR-detected X-ray AGNs in the Chandra Deep fields, split in ter ms of intrinsic rest- 
frame 2-10 keV luminosity (i.e., Lx; see key). Stellar masses are calculated by fitting the optical-NIR fluxes with the stellar templates of lBruzuaT & Charloll 
120031) . We have ensured that the light in these bands is dominated by the host galaxies, rather than the AGN (see i|4.2l and fig- [3}. The distribution of host 
galaxy masses for AGNs of Lx= 10 42 -10 44 ergs s~ 1 is remarkably flat across the considered redshift range. Middle: Same as left panel, but for all host galaxies 
for which we can measure reliable stellar masses (see i]4.2t . Right: Same as left, but now we show the mean galaxy masses in the z = 0.5 — 1, 1—2 and 2 — 3 
bins (separated into Lx bins; see key). Again, there is no clear evidence of any significant evolution in the stellar masses of AGN host galaxies out to z < 3. In 
all plots we also show the distribution of masses for the underlying galaxy population (calculated in the same way as for the AGN hosts) and in the right-hand 
plot we include th e values of M* (i.e., the point at which the galaxy mass function turns over) for star-forming galaxies as a function of z as reported in 
lllbert et al.l JioTol) . 



nosities of Lx= 10 42 -10 44 ergs s _1 AGNs is consistent with the 
increase in infrared output of normal, star-forming (hereafter, main- 
sequence) galaxies. They argue that this could imply that moderate 
luminosity AGNs typically reside in main-sequence galaxies under- 
going internally-induced evolution (as opposed to major-mergers). 
However, there remains other possible explanations for the increase 
in Lir of these AGNs; changes in the average mass of AGN-hosting 
galaxies or a transition from starbursting to main-sequence hosts 
could also produce similar trends. In this section, we will combine 
the results outlined above with mass estimates of the host galaxies 
(see Q4.2\ to discriminate between these three scenarios. 

To determine whether the increase in Lir of moderate lu- 
minosity AGNs is due to an increase in the average host galaxy 
mass towards higher redshifts, we plot the stellar masses of the 
AGNs in our sample as a function of redshift (see fig. |7J. In the 
left-hand panel of fig. [7] we plot those Lx= 10 42 — 10 44 ergs s -1 
AGNs that are detected at FIR wavelengths on top of the under- 
lying distribution of stellar masses of non-AGNs in the GOODS 
fields. We find no significant evolution in the stellar masses of 
FIR-detected AGN hosts over the redshift range considered (i.e., 
0.5 < z < 3), with these moderate luminosity AGNs residing 
almost exclusively in 10 10 — 10 12 M Q galaxies at all redshifts 
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10 44 ergs s~ - 1 AGNs). IXue et alj ( bold) reported similar findings 
for L x ~ 10 42 -10 44 ergs s" 1 AGNs in the CDFs out to z ~ 4 
and this mass range is also similar to that of optical narrow- 
line selected AGNs in the local Universe (i.e., 0.02 < z < 0.3; 
lKauffmannetal.ll2003l) . We get marginally smaller mean stellar 
masses when we consider all AGNs for which we can calculate 
secure host galaxy masses (i.e., ~69 per cent; see right-panel of 
fig.|7J <M sta , s >=8.6+J g x 10 10 Mq), indicating that detection at 
FIR wavelengths biases in favour of slightly higher host galaxy stel- 
lar masses (although we note that the mean masses are consistent 



to within la uncertainties). Interestingly, the AGNs in our sample 
show a weak correlation between Lx and host galaxy stellar mass, 
which can be expressed as: 



star 6xl0 4 



ergs s 



1/7 



(4) 



This corresponds to a change in the median stellar mass of only a 
factor of ~ 2 across two orders of magnitude change in Lx (i.e., 
Lx= 10 42 -10 44 ergs s _1 ; see fig. [8} and is too weak to account for 
the prevalence of AGNs among only high mass galaxies01n other 
words, the high host galaxy stellar masses of our AGN sample is 
not the result of a selection effect. We stress that these results only 
apply to AGNs with Lx= 10 42 -10 44 ergs s _1 , it is entirely possi- 
ble that lower luminosity AGNs (including those below the detec- 
tion threshold of the CDF surveys) reside in lower mass galaxies. 
However, since there is no significant evolution in the host galaxy 
stellar masses of moderate luminosity AGNs we can rule out the 
possibility that the observed increase in their average Lir is due 
to a systematic increase in their host galaxy stellar mass with red- 
shift. Instead, we must attribute it to an increase in their SSFRs. 
Next, we consider whether this increase in their SSFRs is con- 
sistent with that observed in the general star-forming (i.e., main- 
sequence) galaxy population, or whether it signifies a transitioning 
of AGN host galaxies from starbursting galaxies at z = 2 — 3 to 
main-sequence/quiescent galaxies today. 

In the formalism of El baz et alj d20 lTh . secure starburst galax- 
ies have SSFRs that are > 3 times higher than those of main- 
sequence galaxies (here, "secure" implies avoiding any contamina- 



10 We note that this correlation remains largely unchanged when we ex- 
clude AGNs from our highest redshift bin (i.e., z = 2 — 3), suggesting that 
selection biases do not have a major effect. 
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Figure 8. Intrinsic (absorption-corrected) rest-frame 2-10 keV X-ray lu- 
minosity (i.e., Lx) vs. host galaxy stellar mass (i.e., M stMS ) for the X-ray 
AGNs in the Chandra Deep Fields, separated in terms of redshift (see 
key). The dashed lines indicate the Lx= 10 42 -10 44 ergss -1 range used 
throughout this study. Large black stars are used to show the median stel- 
lar masses of AGNs in six different Lx bins (0.5 dex wide, starting at 
log (Lx/ergs s _1 ) = 41.5. The solid line is a least-squares linear fit to these 
six points. We find only a weak correlation between Lx and host stellar 
masses for the AGNs in our sample; over two orders of magnitude increase 
in Lx (i.e., Lx= 10 42 -10 44 ) the median stellar mass only increases by a 
factor of ss2. 

tion from main-sequence galaxies). This is the case at all redshifts, 
although the SSFRs of main-sequence galaxies increase at higher 
redshifts. Thus, to determine whether the observed Ljr-z trend for 
X-ray AGNs is caused by a transitioning of their hosts from star- 
bursting to main-sequence galaxies, we will compare their SSFRs 
over the 0.5 < z < 3 redshift interval of our sample. 

In fig.|9]we plot the SSFRs of the L x = 10 42 -10 44 ergs s" 1 
AGN-hosting galaxies in our sample as a function of redshift. In- 
cluded in these plots are the SSFRs of normal (i.e., non-AGN- 
hosting) galaxies withM stars = 10 10 -10 12 M Q detected in GOODS- 
H, i.e., only those galaxies with stellar masses within the range of 
the AGNs in our sample. We also include tracks sho wing the SSFR- 
z trend s fo r main-sequence g alaxies reported in iPannella et al.l 
J2009ij) and lElbaz et all d20 1 J see their Eq. 13)0 In the left hand 
panel of fig.|9]we plot the SSFRs for the FIR-detected AGNs in our 
sample. The majori ty (i.e., ~81 per ce nt) of these AGNs lie within 
a factor of 3 of the lElbazetaljBoTIt) main-sequence trend, with 
roughly equal numbers of FIR-detected AGNs outside this range 
having higher or lower SSFRs (i.e., ~47 per cent and ~53 per cent 
of outliers, respectively). However, none of the FIR-det ected AGNs 
in our sample have SSFRs lower than 17 per cent of the Elb az et ail 
J201 l|) main-sequence and a re thus unlikely to be quiescent galax- 
ies (c.f.. lNoeskeet"aT] l2007l). Later in this section we use the re- 



11 A fit to the SSFRs of the comp arison M sais = 10 10 -10 12 M Q galaxies 
shown in fig.[9]reproduces Eq. 13 o f lElbazetalJfeOllh . 



suits from our stacking analyses to estimate what fraction of all 
(i.e., FIR detected and undetected) L x = 10 42 -10 44 ergs s _1 AGNs 
reside in quiescent galaxies. Over the entire 0.5 < z < 3 redshift 
range only 8 (i.e., ~10 per cent) of the 77 FIR-detected, Lx= 10 42 - 
10 44 ergs s" 1 AGNs with secure host galaxy stellar masses have 
SSFRs typical of starburst galaxies at their redshifts. This number 
rises to 1 1 when we include the 3a upper limits of FIR-undetected 
sources (i.e., ~ 8 per cent of all 130 AGNs within these redshift 
and Lx bounds). 

The trends seen in the FIR-detected fraction of our AGN sam- 
ple are reflected in the average SSFRs derived from our stacking 
analyses (see fig. [9] right-hand panel). The average S SFRs of the 
AGN host galaxies in our sample lie slightly below the Elb az et all 
d201ll) trend for main-sequence galaxies. At z = 0.5 — 1, 1 — 2 
and 2-3, L x = 10 42 -10 44 AGNs have average SSFRs that are 
71^ 2 g, 89^ 2 g and 72^ 1 P er cen t of those of typical main-sequence 
galaxies at their average redshifts, respectively (i.e., z =0.81, 1.4 
and 2.5)0 

Since we are effectively comparing against the aver- 
age SSFRs of M star s= 10 10 — 10 12 M© galaxies (i.e., see footnote 
lilt , it is unlikely that these lower SSFRs are due to the AGN hosts 
having higher-than-average stellar masses. Despite being slightly 
lower, the mean SSFRs of the AGNs in our sample show largely 
the same magnitude increase as main-sequence galaxies, irrespec- 
tive of Lx; a factor of 3.1^ 2 'J from z =0.81 to 1.4 and a factor of 
7.8±f;^ from z =0.81 to 2.5 forL x = 10 42 -10 44 ergs s _1 AGNs. 

As we have already seen, none of the FIR-detected AGNs in 
our sample have SSFRs typical of quiescent galaxies (i.e., < 10 per 
cent of the average SSFR of main-sequence galaxies at a given red- 
shift). However, this still leaves the ~60 per cent of AGNs in our 
sample that remain undetected at FIR wavelengths. We can use the 
results from our stacking analyses to place constraints on the frac- 
tion of these FIR-undetected, L x = 10 42 -10 44 ergs s -1 AGNs that 
reside in quiescent galaxies. Adopting a Monte Carlo approach, we 
randomly assign a < 3o flux to each FIR-undetected source as- 
suming an underlying log-normal SSFR distribution (Rodighiero 
et al. in prep.), then calculate the average SSFR of these simu- 
lated sources. We repeat this process 1000 times, selecting only 
those trials whose simulated average SSFRs lie within the mini- 
mum and maximum bounds derived from our stack. For each suc- 
cessful trial, we calculate the fraction of AGNs with simulated 
SSFRs below 10 per cent of SSFRms (i- e -> quiescent). In doing so, 
we estimate that between 13 and 37 (i.e., ~16 and ~46 per cent, 
respectively) of the 81 FIR-undetected AGNs between z = 0.5 — 3 
with securely-determined stellar masses have SSFRs that are con- 
sistent with those of quiescent galaxies. These numbers correspond 
to 15 ± 7 per cent of all L x = 10 42 -10 44 ergs s -1 AGNs within this 
redshift interval. This fraction is similar in size to the average deficit 
of AGN SSFRs compared to the main-sequence trend, suggesting 
that the inclusion of quiescent galaxies in our average SSFRs is 
responsible for this offset (see also 36j};EB By incorporating the 
results derived from the FIR-detected AGNs, we estimate that 7±2 
per cent of L x = 10 42 -10 44 AGNs at z =0.5-3 reside in starburst- 
ing galaxies and 79 ± 10 per cent reside in main-sequence galax- 



12 The corresponding results for the Lx= 10 42 -10 43 ergs s 1 bins are 
75^ 4 , 87±|| and 91±l\ per cent; for the L x = lO^-lO 44 ergs s -1 bins they 
are 6Q±H, 98±|| and 73+2? P er cent - 

13 By excluding the quiescent 15±7 per cent of hosts from our average 
SSFRs, we estimate that the SSFRs of star-forming AGN hosts are ~95 
per cent of those of typical main-sequence (i.e., non-AGN) galaxies (and 
consistent to within the various uncertainties). 
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Figure 9. Left: The specific star formation rates (i.e., SSFRs) of the host galaxies of the X-ray AGNs in our sample, plotted against redshift and separated 
in terms of intrinsic rest-frame 2-10 keV luminosity (see key). Downward pointing arrows are used to indicate the 3o upper limits on SSFR in cases where 
there is no formal FIR detection. The vast majority of FIR-detected AGN reside in normal, star-forming (i.e., main-sequence) galaxies. Right: Same as the 
left panel, but for the mean SSFRs calculated using the mean star formation rates (i.e., SFRs) and stellar masses (i.e., M stilrs ) of the AGNs in various redshift 
and Lx bins (see key). Also included in this plot are the average SSFRs calculated from the stacks of the FIR-undetected AGNs which we use to estimate the 
fraction of AGNs that reside in quiescent galaxies (solid black triangles; see i]5.2t . For comparison, we include in both panels the non-AGNs detected in the 
GOODS-H survey (restricted to M sta is= 10 10 — 10 12 Mp. to match the range of stellar masses spanned by our AGN sample; see i]5.2l and fig. [7). Also included 
are tracks showin g the SSFR s of normal, s tar-forming (i.e., main-sequence) galaxies at the redshifts considered (shown by the solid and dashed lines; from 
lElbaz et alj[201 ll and lPannella et alj|2009ij, resp ectively). The do tted line is a factor of 3 above the Elb az et alj 1201 11) curve, above which a galaxy can be 
safely considered to be strongly starbursting (see lElbaz et alj201lh . 



ies. These fractions do not change significantly with redshift and 
show a slight preference toward main-s equence galaxies com pared 
to the general population. For example. iNoeske et alj (12007) found 
that ~56-67 per cent of 0.2 < z < 1.1 galaxies (with similar stel- 
lar masses as the AGN hosts considered here) are main-sequence, 
while ~3 per cent are quies cent and ~1 per cent are starbursts 
(although Noesk eet alj d2007f) defined starbursts as having SSFRs 
>5 times higher than that of main sequence galaxies, rather than 
>3 times used here). 

Whilst we can only place upper limits on the fractions of 
AGNs that reside in quiescent and main-sequence galaxies, we are 
able to report that the SSFRs outlined here are not compatible with 
a view that AGN hosts tend to be strongly starbursting galaxies, 
nor that moderate luminosity AGNs have transitioned from resid- 
ing in starburst galaxies at high redshifts to m ain-sequence gala xies 
at z ~ 0.5. This is in general agreement with lXue et alj J2010h . al- 
though they found that z < 1 X-ray AGNs in the CDFs have SFRs 
that are a factor of 2-3 times higher compared to a stellar mass- 
matched sample of non-AGNs. However, we note that their SFRs 
were derived from 24 |j.m observations, which are more suscepti- 
ble to contamination from either emission from the AGN or strong 
spectral features (i.e., PAH lines). 

Our results disfavour changes in either the stellar masses 
of AGN hosts and/or a transitioning of their hosts from main- 
sequence/quiescent to starbursting galaxies as giving rise to the 
observed increase in AGN Lir at higher redshifts. We have also 
shown that the SSFRs of AGN hosts are only marginally lower 



(i.e., ~ 20 per cent lower), than those of main-sequence galaxies at 
0.5 < z < 3. This extends the work of lLutz et alj (2010) by quanti- 
fying their suggestion that moderate luminosity AGNs typically re- 
side in main-sequence galaxies undergoing internally-induced evo- 
lution. We further explore the consequences of this interpretation 
inU 

5.2.1 The red rest-frame colours of AGN -hosting galaxies 

Here we take a brief aside to address a specific question regarding 
the rest-frame optical colours of the AGN host galaxies explored 
in this study: if the majority of moderate-luminosity AGNs reside 
in star-forming, main-sequence galaxies, then why are their host 
galaxies generally redder than usually expected for these types of 
galaxies (i.e, see fig-O? The ~ 20 per cent difference in the average 
SSFRs of AGN hosts is too small to account for their red colours. 
Instead, the answer likely lies in the high masses of AGN host 
galaxies. More massive star-forming galaxies tend to be redder than 
their lower mass counterparts with similar SSFRs d ue to increased 
levels of dust obscuration in these galaxies (e.g., Pan nella et alj 
2009a). Indeed, when we compare the rest-frame colours of our 
(L x = 10 42 -10 44 ergs s _1 ) AGN hosts against those of non-AGN- 
hosting, main-sequence galaxies with similar stellar masses (i.e., 
within ±0.2 dex of the median stellar mas s of the AGN ho sts), we 
find that they are remarkably similar (see IXue et af]|201dl and our 
fig. |10t. This is contrary to a number of recent studies reporting that 
X-ray AGNs often have rest-frame colours offset from the majority 
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Figure 10. Same as fig. [5] but now we compare the colours of our sample of X-ray AGNs against only Herschel-detected main-sequence galaxies with stellar 
masses (i.e., M stals ) within ±0.2 dex of the characteristic M sta rs of AGN host galaxies (i.e., 1 .53 x 10 1 1 Mq). The host galaxy colours of moderate luminosity 
(i.e., Z-x= 10 42 -10 44 ergs s _1 AGNs is remarkably similar to those of main-sequence galaxies with similar stellar masses (see £|5.2.U . 
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Figure 11. Plot of rest frame U — V vs. host galaxy stellar mass for the 
AGNs in our sample (large symbols) and non-AGNs in the same fields 
(small symbols). The rest frame colours have not been corrected for the 
effects of reddening. There is a tendency for the AGNs in our sample tend 
to reside in either the so-called "red sequence" or "green valley", despite 
their far-infrared derived SFRs and SSFRs indicating that the majority are, 
in fact, star-forming (i.e., main-sequence) galaxies. We conclude that the 
reddened colours of these massive, star-forming AGN hosts are largely due 
to the effects of dust obscuration. Here, we have included all Herschel de- 
tected AGNs, not just those with "secure" masses (see i|4.2l . so as not to 
bias the plot. 



of the galaxy population (i.e., within the so-called "green valley", 
rather than either the "red sequence" or "blue cloud"; see our fig 
lilt , which has been interpret ed as evidence of the quenching of star 
formation by the AGN (e.g. Nandra et alj 2007 : Schawins ki et al 



20071: iGeorgakakis et alJl2008l: iBundv et alfeooa l Silverman et al 
20081) . Instead, once the effects of reddening are correctly taken 
into account, AGN hosts display largely the same distribution of 



intrinsic colours as the general ga laxy population, as was recently 
shown bv lCardamone et aLH20ld) . As further confirmation of this, 
the fractions of main-sequence and quiescent AGN hosts reported 
here (see t]5.2| l are broadly consiste nt with the fractions of star- 
forming (blue cloud) hosts reported in ICardamone eTail d2010h . al- 
though they report slightly higher fractions of quiescent (red se- 
quence) hosts (likely due to the different approaches used to iden- 
tify quiescent galaxies, i.e., infrared vs. optical diagnostics). Thus, 
the hosts of moderate AGNs have largely the same fundamental 
stellar properties (i.e., rest-frame colours, SFRs, SSFRs) as their 
non-AGN-hosting counterparts of similar stellar mass. 



5.3 The optical to far-infrared SEDs of X-ray AGNs 

If the observed increase in Ljr of AGNs with redshift is solely due 
to the evolving SSFR of their main-sequence hosts, as our results 
seem to suggest, we should expect their infrared SEDs to be in- 
creasingly dominated by star formation at higher redshifts. In fig. 
[T2I we plot the average rest-frame SEDs traced out by the opti- 
cal to far-infrared bands for Lx= 10 42 — 10 44 ergs s _1 AGNs in 
each of our redshift bins. While the optical to near infrared por- 
tions of the average AGN SEDs show only modest change with 
redshift, both the mid (i.e., 16 and 24 p.m) and far-infrared (i.e., 
100 and 160 portions show a strong increase with i ncreasing 
redshift, as previously reported by Mullane y et alj d2010T) . This in- 
crease is strongest at FIR wavelengths meaning that the average 
infrared SED of L x = 10 42 -10 44 ergs s~' AGNs is redder at higher 
redshifts. 

Is the change in the average AGN SED consistent with pure 
evolution of the SSFR of main-sequence galaxies? To assess this 
we predict the infrared SEDs of high redshift AGNs by crudely 
assuming that only the SSFRs of their host galaxies have evolved 
with redshift, then compare the observed SEDs against these pre- 
dictions. To generate the predicted SEDs we first fit the average in- 
frared SED of z ~ BAT/IRAS AGNs with a model comprising an 
AGN and a single host galaxy component (see Mullane y et alj201ll 
for a description of the fitting procedure). The AGN component de- 
rived from this fit is then kept constant throughout. To generate the 
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Figure 12. Average optical to far-infrared SEDs of the Lx= 10 42 — 
10 44 ergs s" 1 AGNs in our three moderate/high redshift bins (see key). The 
strong increase in the average infrared luminosities of the AGNs is clearly 
evident in these SEDs. A strong stellar component is evident in all three 
SEDs around 1 |xm, demonstrating that, on average, the host galaxy dom- 
inates the optical and near-infrared fluxes. As further confirmation of this, 
there is even evidence of the 1.6 |J.m stellar "bump" in the average SED 
of the z = 2 — 3 AGNs (where the mean redshift of the AGNs places this 
feature within the IRAC 5.6 |xm filter). 



host-galaxy components for the high redshift AGNs, we estimate 
the expected Lir of the host s based on their averag e stellar masses 
and redshifts (usi ng Eq. 2 of Pannella et "ai]|20 Q9a). We then select 
the SED from the lCharv & Elbaj ( 1200 lh library that most closely 
matches this Lir and add it to the AGN component from the z ~ 
fit to produce the total predicted SED (i.e., AGN+host). The SEDs 
generated in this way are shown in fig. [13] and are roughly consis- 
tent with the observed mean AGN SEDs in each redshift bin (the 
slight discrepancy at MIR wavelengths for our z = 2 — 3 bin is likely 
due to the increased strength of PAH features in high -z star-forming 
galaxies compared to th e local templates used here: ICharv & Popa 
l201cilElbazetal.ll201lh . Since these predictions only take into ac- 
count the change in the SSFR from z ~ 0, we may not expect to see 
any evidence of the reduced SSFRs of AGNs compared to main- 
sequence AGNs (see i]5.1.2t in this plot, especially if z ~ AGNs 
also have similarly reduced SSFRs (e.g.. lSalim et alj2007l) . Whilst 
this is only a crude analysis, it does show that the whole AGN+host 
galaxy infrared SED, rather than just the FIR portion, evolves in a 
manner broadly consistent with an evolving main-sequence SSFR. 
An important prediction that arises from this scenario is that the 
average MIR spectra of AGNs are increasingly dominated by host 
galaxy spectral features (e.g., PAH features) at higher redshifts. We 
will use archival Spitzer-IRS spectra to determine whether this is 
the case in a future paper. 



6 DISCUSSION 

We have used the deepest FIR data currently available to probe 
the SFRs and SSFRs of AGN host galaxies to z < 3. Our main 
results are: (a) there is no correlation between star formation and 
Lx in moderate luminosity AGNs (i.e., Lx= 10 42 — 10 44 ergs S _1 ), 
(b) there is a strong, monotonic, increase in the SFRs and SSFRs 



of the hosts of these AGNs from z ~ to 3 (c) that < 10 per cent of 
AGNs reside in starbursting galaxies, the majority being hosted by 
main-sequence, star-forming, galaxies and (d) the increase in the 
SSFRs of the AGN hosts is broadly consistent with the increase in 
the SSFRs of main-sequence galaxies in general. In this section, 
we explore the possible implication of these results in terms of our 
understanding of the co-evolution of AGNs and their host galaxies. 



6.1 The influence of moderate luminosity AGNs on star 
formation 

The lack of any strong correlation between Lx and Lir in any of 
our redshift bins suggests that instantaneous nuclear and global star 
formation activity are largely decoupled for moderate luminosity 
AGNs at z < 3. At first sight, this would seem to imply that nuclear 
activity has little impact on the levels of star formation taking place 
within their host galaxies, bringing the role of AGN/galaxy feed- 
back in the evolution of these systems into question. However, it 
has recently been reported that nuclear (i.e., < 1 kpc), rather than 
global, SFRs are strongly co rrelated with black hole accretio n rate 
and, thus, AGN luminosity (Dia mond-Stanic & R ieke 2011). This 
would suggest that the physical links between black hole growth 
and star-formation activity are localised to central regions of the 
host galaxy. We should also consider the effects that X-ray vari- 
ability will have on our results. Since X-ray observations only take 
a snapshot of nuclear activity, the lack of a correlation between Lx 
and Lir could simply be due to AGN variability introducing scat- 
ter (e.g.. lHopkins & Hernquistll2006l) . One way to test this, which 
is beyond the scope of this paper, would be to use X-ray stacking 
analyses to calculate the average Lx of all galaxies as a function 
of Lir, assuming that a wide cross-section of the galaxy population 
will sample all stages of AGN variability at a given epoch. 

The evidence that the SSFRs of galaxies hosting moderate 
luminosity AGNs are marginally (i.e., ~ 20 per cent) lower than 
their non-AGN-hosting counterparts, whilst tentative, may hint to 
a causal connection between nuclear activity and star formation in 
their hosts. However, with the data used here, it is difficult to deter- 
mine whether the AGN is playing a direct role in the transitioning 
of a galaxy from the main-sequence to a more quiescent state (i.e., 
quenching). As we suggested in S |5.2| it may simply be the case 
that a non-negligible fraction of moderate luminosity AGNs are 
triggered in quiescent galaxies, rather than in star-forming galaxies 
galaxies which then transition to a more quiescent state. Besides, 
the difference is so slight (and in some of our redshift bins, neg- 
ligible) it is difficult to argue that quenching by moderate levels 
of nuclear activity has any significant impact on the star-formation 
activity of their host galaxies. 

As most previous studies of AGN hosts have largely focussed 
on SFRs, rather than SSFRs, it is difficult to directly compare our 
findings with these earlier studies. In any case, th ey have arrived 
at conflicting con clusions, reporting increased (e.g.jGu et alj200rj . 
lLacv et alj |2007). suppressed (e.g.. lHoll2005t iKim et al. 2006) or 
no significant difference (e.g.. ISilverman et alj|2009f) in the SFRs 
of AGN hosts depending on the type/luminosity of the AGNs 
in their samples and/or selection method. By focussing on SFRs 
(rather than SSFRs) these earlier studies were susceptible to differ- 
ences between the masses of their AGN and control (galaxy) sam- 
ples. However, by comparing the SFRs of mass- matched samples 
of galaxies (synonymous to comparing SSFRs). [Xue et aL IHouJ) 
found that the SFRs of AGN hosts at z = - 1 are a factor of 2-3 
higher than those non-AGNs, whereas this difference largely dis- 
appears at z = 1 — 3. By contrast, in their study of the ultraviolet- 
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Figure 13. The average observed infrared spectral energy distributions (SEDs) of the AGNs in the Chandra Deep Fields and the local BAT//SAS sample 
(points). Again, we have separated our sample in terms of X-ray luminosity (left and right panels; see labels). Also shown on these plots is the expected 
SEDs calculated by boosting the host galaxy component of the fit to the z ~ observed SED by a factor equivalent to the increase in the SSFR at the average 
redshift of each bin (solid lines, colour corresponding to the appropriate redshift bin; see i|5.3l for a full description). The dashed black line indicates the AGN 
co mponent at all re dshifts, while the dotted lines (various colours) represent the host galaxy component at the mean redshift of our subsamples (taken from 
the Chary & Elbaz 2001 library). It is important to note that the solid lines are not an attempted fit to the data, but the construct of a model based on the simple 
assumption that only the normalisation of the host-galaxy component increases with redshift. This plot demonstrates that the increase in luminosity observed 
in each infrared band (i.e., 16 |xm, 24 |xm, 100 |xm and 160 |xm) is roughly consistent with the increase in SSFR observed in normal, star-forming galaxies 
(see 631 . 



derived SSFRs of galaxies in the local Universe, ISalim etalj ( l2007h 
found that AGNs have suppressed SSFRs compared to star-forming 
galaxies, although they found that AGNs had, on average, ~ 60 per 
cent lower SSFRs, compared to the ~ 20 per cent reported here 
(see their Figs. 17 and 18). At face value, this latter point may sug- 
gest that nuclear activity plays a stronger role in quenching star- 
formation at low redshifts. However, it is important to note that 
our sample selecti on and analyses differ considerably from that of 
ISalim et al.l ( 120071) . which may account for some of the disagree- 
ments between these studies. In the end, clear evidence of AGN - 
driven feedback (e.g-. lNesvadba et alj2006j ; I Alexander et alj2010l) 
will be needed to demonstrate that any suppression of star forma- 
tion is due to them having a direct, causal influence on their host 
galaxies. 



6.2 The majority of moderate luminosity, distant AGNs 
reside in normal, star-forming galaxies at all redshifts. 

Our results provide convincing evidence that the majority of dis- 
tant, moderate luminosity AGNs reside in main-sequence, star- 
forming galaxies. To interpret the implications of this result in 
terms of the co-evolution of AGNs and their host galaxies, we will 
consider the typical characteristics of galaxies that lie both on and 
off the main-sequence. 

Main-sequence galaxies tend to be late-type spiral s undergo- 
ing internally driven evolution (e.g. jNoeske et al]|2007h . Star for- 
mation in these galaxies predominantly takes place in their disks 
and is likely triggered by internal processes such as turbulence an d 
disk instabilities (e.g.. [Elmegreen et alj|2007l ; iGenzel et alj|2008h . 



In contrast, starbursting galaxies - the other main category of star- 
forming galaxies - often show signs of having recently under- 
gone a major merger event (mass ratio < 3 : 1), meani ng their star 
formation is more likel y to be merger induced (e.g., iDaddi et all 
120071 : lElbaz et al.ll2007h . It is important to note that the starburst 
phase of a typical major merger is relatively short lived, with 
only ar ound V3 of the life of the merger being spent in this 
regime (Di Matte o et al . 2008). However, for the rest of the time, 
the merger is not actively accelerating star formation beyond what 
would be expected under normal, internal processes. Finally, at the 
other extreme are quiescent galaxies that have ceased to be ac- 
tively for ming stars and have SSFRs significantly lower (i.e., <C 10 
per cent; iNoeske et alj|2007h than those of typical main-sequence 
galaxies. 

Based on this emerging picture, the results presented here 
strongly suggest that moderate levels of nuclear activity have 
mostly taken place in galaxies undergoing internally induced, low 
efficiency star formation, as opposed to starbursting m ergers - a 
view consistent with the conclusions of lXue et al.l (2010). Of course 
a fraction of these AGNs may reside in merging systems that have 
ceased starbursting, although even a conservative estimate places 
this fraction at around 21 ± 6 per cent0 However, even if this is 
the case, any merging has since ceased triggering starbursts, mean- 
ing the merger event is no longer having a significant effect on 
the star formation activity of these galaxies. Furthermore, since it 

14 Based on 7 ± 2 per cent of the AGNs in our sample having starburst-level 
SSFRs, and assuming that a typical merger spends only 1 of its time in a 
starburst phase. 
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seems unlikely that a merger can channel material onto a black hole 
without also inducing significant amounts of star formation, this 
interpretation also implies that the majority of moderate nuclear 
activity has been triggered and driven by internal processes since 
z ~ 3. Since main-sequence galaxies at high redshifts (i.e., z ^ 1) 
are often dominated by gravitationally unstable, clumpy disks (e.g. , 
Cowie et al]|l996l: lElmegreen et af] 12004 lElmegreen et al J 120071 : 
Fdrster Schreiber et alj|200d 120091 ; iGenzel et alj|2008l) . we spec- 



ulate that disk instabilities are driving the inflows that lead to nu- 
clear activity in their cores. Indeed, theoretical models have show 
that the giant clumps can directly fuel black hole growth (e.g., 
lElmegreen e"t~aU 2008), and can more generally trigger a central 
gas inflow fuelling an AGN (Bournaud & Dekel, in prep.). 

These findings are also consistent with the conclusions of 
recent studies that have looked at recent merger histories and 
morphologies of galaxies hosting moderate luminosity AGNs. 
For example, when they compared the morphologies of z < 1 
AGN- hosting galaxies with those of the general galaxy popu- 
lation, Icisternas et al.l d201ll) showed that the same fraction of 
each sample showed si gns of recent or ongoing mergers, while 
ISchaw inski et al] d201ll) used a subsample of 28 X-ray AGNs in 
the CDFS to show that this was also likely the case at z < 3 (see 
also Kocevski et al, submitted, who have performed morphologi- 
cal analyses for a larger sample of AGNs, many of which overlap 
the sample used in this study). Our results build on these earlier 
findings by using infrared-derived SSFRs - arguably a more reli- 
able and quantifiable means of identifying main-sequence galaxies 
- to show that the vast majority of all moderate nuclear activity in 
these deep fields is accompanied by internally-induced, rather than 
violent merger-induced star-formation. 

Whilst this interpretation appears to apply for moderate lumi- 
nosity AGNs at all redshifts since z ~ 3, it is particularly pertinent 
at z < 1, when these AGNs dominated the integrated black hole 
growth taking place in the Universe. Indeed, integrati ng the ana- 
lytic AGN X-ray luminosity functions of lUeda et al.l ( 120031) and 
assuming the lumino sity-dependent bolometric correction factor of 
iMarconi et ail d2004 reveals that L x = 10 42 -10 44 ergs s" 1 AGNs 
have accounted for ~65 per cent of all black hole growth since 
Z ~ 1. Therefore, it would seem that the majority of the volume- 
averaged black hole mass accretion since z < 1 has taken place in 
galaxies undergoing internally-induced evolution, which also ap- 
pears to be the case for the majority of stellar mass build u p over 
this redshift ranee (e.g.. lBell et aHl2005l:|jog~ee et all 2009). What 
this means for the concurrent build up of black hole and stellar 
mass to today's MsH~^Bu\ze relation will be explored in a follow- 
up paper (Mullaney et al., in prep.) but it is important to bear in 
mind that, by assuming the same luminosity functions and bolo- 
metric corrections (although see footnote |2), we calculate that the 
majority (i.e., ~ 60 per cent) of all black hole mass accretion (i.e., 
growth) in the Universe since z ~ 3 would have taken place in more 
luminous AGNs (i.e., Lx^ 10 44 ergs s -1 ). 

On the other hand, evidence supporting the importance of 
mergers on the concurrent growth of black holes and their host 
galaxies has been present e d in a number of previous studies. For 
example, Alexa nder et al.l «j2005) showed that ~ 30 — 50 per cent 
of sub-mm galaxies, which often lie above the main-sequence 
(in terms of SSFR), contain moderate to high luminosity (i.e. , 
log[L x /ergs s~'l > 43.5) AGNs (see also, e.g., iLutz ,et alj|2005t 
Valian te et al J 12007]: iMenendez-Delmestre et alj|2007l: Pope et all 
20081 : ICoppin et aljboiOt) . Furthermore. IPage et al.l d2004 showed 
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that a significant proportion of luminous, absorbed quasars reside 
in strongly star-forming galaxies (although some of these may now 



Figure 14. Top: The space density of moderate luminosity AGNs as a func- 
tion of host galaxy stellar mass (i.e., M slars ). These space densi ties ar e cal- 
culated by multiplying the AG N fractions report ed in lXue et al . (2010) with 
the total mass functions from Ilbert et alj d2010h . The m ass functions of all , 
quiescent and star-forming galaxies are also taken from Ilbert et al. ( 2010; 
see key). Both AGN mass functions peak at roughly the point where the 
total galaxy mass function turns over (i.e., M t ). Bottom: The fractions of 
main-sequence, star-forming galaxies that host moderate luminosity AGNs 
as a function of stellar mass. The fraction of main-sequence galaxies host- 
ing moderate luminosity AGNs increases strongly with host stellar mass. 
This, in turn, implies that the AGN duty cycle also increases strongly with 
host galaxy stellar mass. 



be recog nised as main- s equen ce galaxies at high redshifts). More 
recently, Ferugl kTet al.l d201ll) identified two Compton-thick, lu- 
minous (i.e., L x > 10 44 ergs s" 1 ) AGNs in the CDF-S whose host 
galaxies are also undergoing a period of intense star formation. 
Taken in conjunction with the arguments presented above it is 
tempting to speculate that, while internal processes have largely 
driven moderate nuclear activity, external (i.e., non-secular, major- 
merger driven) processes have played a more major, perhaps dom- 
inant, role in the build up of black hole and galaxy bulge mass 
and, thus, formed today's M^-M^ ^ relationship - as su ggested 
by the theoretical models o f Hopkins & Hernquisj d200fj|) and re- 
cently discussed in ILutz et al.l fcOlCh . However, this interpretation 
would be difficult to reconcile with emerging evidence that long 
lived, internal processes were the main driver for st ar formation at 
these high red s hifts ( e.g. jDaddi et alJl2008Ll2010ah . Furthermore, 
iMainieri et al.l d201lh recently used a combination of optical spec- 
troscopy and SED fitting to measure the host-galaxy properties of 
luminous (i.e., L x > 10 44 ergs s _1 ) Type-2 quasars, reporting re- 
markably similar fractions of main-sequence and quiescent hosts 
as those found here for moderate luminosity AGNs. If confirmed, 
such findings would imply that the majority of all nuclear activity 
has taken place in galaxies undergoing internally-driven evolution 
and not in major-mergers. 
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6.3 The fraction of main-sequence galaxies that host AGNs 

We have shown that the majority (i.e., 79 ± 10 per cent) of dis- 
tant, moderate luminosity AGNs reside in main-sequence galaxies. 
However, these AGNs are not evenly distributed among the main- 
sequence population. Instead, they preferentially reside in massive 
galax ies (i.e., M stars > 10 10 M©) b oth locally (e.g., Kauffman n et al 



2003 ) and at high redshifts (e.g., [Silverman et al. 2009; X ue et al 



201C ; also our fig. Therefore, while nuclear activity is rela- 
tively rare in the overall galaxy population (~ 6 — 7 per cent of all 
^stars> 3 x 10 9 M© galaxies), a significant fraction of high mass 
galaxies host moderate luminosity AGNs. We know that locally, up 
to 80 per cent o f M stars = 10 11 M Q galaxies are hosting an AGN 
feest et alj200l) . 

Recently, a number of studies have shown that the mass func- 
tion of main-sequence galaxies has maintained a roughly con- 
stant shape si nce z ~ 2, with only its normalis ation changing 
with redshift dllbert et alJboid : IPeng et alj|20ich . The range of 
characteristic masses at which this mass function turns over (i.e., 
log[M*/M©] = 10.7 — 11.0) is remarkably similar to the typical 
stellar masses of galaxies hosting moderate luminosity AGNs (see 
fig- 12} - Whilst the similarity between these two mass ranges may 
be coincidental, it is interesting to explore it further by considering 
roughly what fraction of main-sequence galaxies within this mass 
range host moderate luminosity AGNs. 

In their recent study of the rest-frame colours, magnitudes 
and host g alaxy masses of AGNs in the Chandra Deep Fields, 
IXue et all feoid) noted that the fraction of galaxies containing 
AGNs in similar Lx and redshift bins to those considered here 
(i.e., z = 0-1, 1-2, 2-3 and log[Lx/ergs s ] =41.9-43.7) increases 
strongly as a function of host galaxy mass, but shows only a weak 
change with redshift. By combining the AGN fractions calculated 
in that study wit h the galaxy stellar mass functions published in 
lllbert et alj i |20ld) . we calculate both the mass function and frac- 
tion of galaxies hosting moderate luminosity AGNs as a function of 
galaxy stellar mass (see fig.U4bP^IWe find that the fraction of main- 
sequence galaxies hosting log[Lx/ergs s _1 ] =41.9-43.7 AGNs in- 
creases strongly from ~l-2 per cent at M stars ~3x 10 9 M© to >20 
per cent at M stars ~3x 10 11 M e , although we note the large uncer- 
tainties in these fractions at high stellar masses due to the rela- 
tively small num ber of AGNs ho sted in M stars > 10 11 M e galaxies 
in the CDFs (see lXue et alfeofol . Fig 12k). Assuming that nuclear 
activity is a largely stochastic process, these fractions imply that 
the duty cycle of AGNs in more massive galaxies is significantly 
longer than in lower mass galaxies. We can rule out selection ef- 
fects as producing this strong increase because (a) there is only a 
very weak correlation betwe en Lx and host stellar mass (see Fig.[8} 
and (b) the lXue e t al. (2010) sample is largely complete within the 
log[L x /ergs s" 1 ] =41.9-43.7 range. Similarly, incompleteness due 
to the effects of X-ray absorption is also unlikely to produce this 
effect. We again stress that these calculations only apply to AGNs 
within the log[Lx/ergs s _1 ] =41.9-43.7 luminosity range. 

Whilst we acknowledge that these fractions are only crude es- 



15 It was recently suggested that stellar mass, combined with emission 
lines, can be used as an effective AGN diagnostic out to z ~ 1 (in the "Mass- 
Excitation" diagram; Juneau et al. 201 1). The lack of evolution in the stellar 
mass of moderate luminosity X-ray AGN since z. ks 3 suggests that the di- 
agnostic power of this parameter could be extended to higher redshift. 

16 We use the m ass function of "intermediate activity" galaxies from 
lllbert et aT] (2010) as they have roughly the same SSFRs as main-sequence 
galaxies at the redshifts considered. 



timates, it is clear that a considerable fraction of main-sequence 
galaxies with masses above the turnover of the galaxy mass func- 
tion host moderate luminosity AGNs. Incidentally, this turnover is 
similar to the stellar mass range that has dominated the star for- 
mation rate density since z ~ 3. While we are reluctant to over- 
interpret this result, especially in light of the arguments outlined 
above, if nuclear activity is necessary to quench the growth of mas- 
sive galaxies, then there is no lack of it in these galaxies. How- 
ever, without clear evidence of AGN "feedback" from these mod- 
erate luminosity AGNs (e.g., powerful, l arge scale outflows suc h as 
those reported in|Nesvadba et al.1120061 ; lAlexander et alj|20ld and 
Rupke & Veilleuxll201lh . we are disinclined to take any firm con- 
clusions on AGN quenching from this analysis. 



7 SUMMARY 

We have used the deepest surveys yet undertaken with the Chandra, 
Spitzer and Herschel telescopes to explore the infrared properties 
and star formation rates of the host galaxies of X-ray AGNs out 
to z ~ 3. By using the ultra-deep GOODS-Herschel observations 
to identify infrared counterparts to X-ray AGNs in the Chandra 
Deep Fields we reach an infrared detection rate roughly 2 times 
higher than any previous far-infrared (i.e., FIR) study of AGNs in 
these fields. In at least 94 per cent of cases of AGNs at z > 0.5, the 
observed 100 |j.m and 160 |J.m wavelengths probed by Herschel- 
PACS are dominated by emission from the host galaxy meaning 
that we can reliably use Herschel fluxes to measure the (specific) 
star formation rates of the host galaxies of AGNs. 

The main results from this study can be summarised as fol- 
lows: 

• There is no clear correlation between AGN power and infrared 
luminosity (and, by inference, star formation rates) for moderate 
luminosity AGNs at all redshifts considered (see ij5.lt . 

• The previously observed increase in the mid to far-infrared lu- 
minosities with redshift is entirely consistent with being caused by 
the increase in the specific star formation rates (i.e., SSFRs) ob- 
served in normal, star-forming (i.e., main- sequence) galaxies (see 

• The majority (i.e 79 ± 10 per cent) of moderate luminos- 
ity (i.e., Lx= 10 42 — 10 44 ergs s _1 ) AGNs reside in main-sequence 
galaxies. We estimate that only 7 ± 1 per cent reside in starburst- 
ing galaxies and 15 ±7 per cent reside in quiescent galaxies. These 
fractions display a slight preference toward main-sequence com- 
pared to the general galaxy population and a remarkably similar to 
those recently reported for Type 2 Quasars (see i]5.2l and §6.2\ . 

• The average SSFRs of moderate AGN hosts are only 
marginally lower (i.e., ~ 20 per cent lower) than those of normal 
main-sequence galaxies, with this small deficit being due to a frac- 
tion of AGNs residing in quiescent galaxies. While this may hint 
at a causal connection between moderate levels of nuclear activity 
and star formation (e.g., AGN quenching ongoing star formation), 
the impact on SSFRs is clearly very limited; (see ^5.2l and [6.1b . 

• The observed 16-160 |i.m SEDs are also broadly consistent 
with being more strongly dominated by the host galaxy at higher 
redshifts (see S j5.3b . Parenthetically, we note that this means that 
using infrared diagnostics to identify AGNs will be increasingly 
difficult at higher redshift. 

• The host galaxies of moderate luminosity AGNs have the 
same red rest-frame colours as non-AGN-hosting galaxies of sim- 
ilar mass. Despite actively forming stars, these massive (i.e., > 
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10 Mq) galaxies have red colours, likely due to their high con- 
centrations of dust (see i ]5.2.1l l. 

• The fraction of main-sequence galaxies experiencing moder- 
ate levels of nuclear activity increases strongly with host galaxy 
stellar mass, from a few per cent at Af s tars~ 3 x 10 9 M© to > 50 per 
cent at M stars ~ 3 x 10 11 M Q (see $63[ . 

In our discussion we argue that the results presented here 
strongly suggest that the majority of moderate nuclear activity in 
the Universe has taken place in normal star forming galaxies under- 
going internal evolution, rather than in violent mergers (see %.2i . 
Since merger-driven inflows are ruled out in the majority of cases, 
we speculate that disk instabilities play an important role in fuelling 
nuclear activity, as suggested by some theoretical models. What is 
remarkable is that we do not see any correlation between the lev- 
els of AGN and star formation activity taking place within these 
galaxies. Indeed, the only significantly different aspect of AGN 
host galaxies compared to the general galaxy population are their 
high stellar masses, which appears to have been the main distin- 
guishing feature of galaxies hosting moderate AGNs since z ~ 3. 
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Table 1. The average AGN and host galaxy properties of X-ray sources in the Chandra Deep Fields, separated into the various classification, Lx and redshift 
bins used throughout this study. 
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NOTES:(l) Index, (2) Description of bin; "All" includes all X-ray detected sources, irrespective of classification (i.e., including AGNs, galaxies and stars), 
rows 4-16 only include AGNs (3) Number of X-ray AGNs with known redshifts covered by the GOODS-H survey, excluding those that show signs of being 

AGN-dominated at FIR-wavelengths (see i]4.U (4) As for column 3, but only counting those that are undetected at FIR wavelengths (i.e., the number of 
sources in each stack). The numbers in brackets are the fractions Mjnd/TVAll (in per cent) (5) As for column 3, but only counting those for which we obtain a 

"secure" host galaxy stellar masses (see i|4.2t . The numbers in brackets are the fractions /VMass /VVaii (in per cent) (6) Mean redshift (7) Mean rest-frame, 
absorption corrected 2-10 keV luminosity (8) Mean stellar mass (9) Mean infrared luminosity calculated from the stacked 100 |xm or 160 (im (depending on 
redshift, see i|4.U of the undetected sources combined with the mean 100 u.m or 160 |xm flux of the detected sources (10) Host galaxy star formation rates 
derived from the infrared luminosities in column 10 using the prescription of Kennicutt ( 1998) (11) Specific star formation rates calculated by dividing the 

star formation rates in column 1 1 by the stellar masses in column 9. 
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